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ABSTRACT 


The  Naval  Postgraduate  School  is  currently  in  possession 
of  software  designed  to  perform  a  thermal  analysis  of 
electronic  components.  This  software  package  incorporates  a 
model  builder  which  contains  two  programs  whose  primary 
function  is  to  generate  a  thermal  model.  In  its  present 
configuration,  the  model  builder  requires  an  inordinate  amount 
of  time  for  data  input  and  model  verification.  This  thesis 
describes  the  development  of  a  model  builder  designed 
specifically  to  reduce  the  time  required  to  model  the 
substrate , epoxy  and  carrier  layers  of  a  microcircuit  assembly. 
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I. 


INTRODUCTION 


The  Naval  Postgraduate  School  is  currently  in  possession 
of  software  designed  to  perforin  a  thermal  analysis  of 
electronic  components.  This  software  package  incorporates  a 
model  builder  which  contains  two  programs  whose  primary 
function  is  to  generate  a  thermal  model,  or  input  data  file, 
to  be  read  by  the  thermal  analyzer  program.  The  first  program 
is  considered  to  be  a  general  model  builder  which  is  used  in 
all  model  development  stages  as  well  as  to  modify  an  existing 
model.  The  second  alternative  was  developed  to  generate  a 
thermal  model  of  a  specific  microcircuit  geometry. 

The  development  of  an  accurate  thermal  model  of  an 
electrical  component  requires  that  the  structure  be  subdivided 
into  a  large  number  of  small  but  finite  subvolumes.  Each 
subvolume  is  assumed  to  be  isothermal  with  the  centroids,  also 
called  nodes,  considered  to  be  r^p-esentative  of  the  entire 
subvolume.  The  most  difficult  problem  encountered  in  the 
development  of  a  thermal  model  is  the  generation  of  «-node 
equations  in  ^-unknown  temperatures  where  the  nodes  are 
connected  by  thermal  conductances.  As  the  desired  accuracy  of 
the  thermal  model  increases,  the  number  of  required  node 
equations  becomes  extremely  large.  Therefore,  it  is  imperative 
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that  the  design  engineer  have  access  to  a  model  builder  that 
will  produce  the  thermal  model  in  a  reasonable  period  of  time. 

In  its  present  configuration,  the  thermal  analysis 
software  contains  a  model  builder  that  generates  the  required 
node  equations  automatically.  There  is  no  question  that  the 
existing  model  builder  programs  have  replaced  the  extremely 
laborious  and  time  consuming  process  of  generating  the  node 
equations  by  hand.  However,  they  still  require  an  inordinate 
amount  of  time  for  data  input  and  model  verification. 

This  thesis  describes  the  development  of  a  model  builder 
designed  specifically  to  reduce  the  time  required  to  model  the 
substrate,  epoxy,  and  carrier  layers  of  a  microcircuit 
assembly.  A  typical  microcircuit  package  configuration  is 
shown  in  Figure  1.  Figure  2a  provides  a  horizontal  interior 
illustration  while  Figure  2b  displays  the  specific  geometry  to 
be  modeled.  All  three  layers  may  contain  an  equal  number  of 
nodes  over  their  width.  However,  the  carrier  layer  may  contain 
a  mounting  ear  on  the  front  and  rear  surfaces.  Additional 
characteristics  to  be  discussed  in  what  follows  are: 

1)  The  capability  of  working  in  English  or  SI  units. 

2)  The  choice  of  four  aspect  ratios. 

3)  The  provision  for  up  to  740  nodes  depending  on  the 

existence  of  mounting  surfaces  (ears)  . 

4)  The  ability  to  input  heat  dissipation  using  several 

methods . 

5)  The  provision  for  six  ambient  temperatures. 
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6)  The  provision  for  rapid,  menu-driven  data  input. 

7)  The  automatic  calculation  of  conductance  values  based  on 
user  input. 


Figure  1.  TO-5  configuration.  (Courtesy  of  Honeywell , Inc . ) 
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II. 


THE  REASON  FOR  THERMAL  ANALYSIS 


Over  the  past  several  decades  a  trend  of  increasing 
sophistication  and  complexity  has  enveloped  the  electronics 
industry.  This  continuing  advance  in  technology  has  greatly 
increased  the  reliability,  capability,  performance,  and 
availability  of  electronic  systems.  The  escalating  demand  for 
further  advances  in  all  areas  of  electronics  has  presented 
engineers  with  an  abundance  of  complex  problems. 

One  major  area  of  concern  is  the  continued  development  of 
advanced  methods  in  the  thermal  control  of  multilayered 
structures.  It  is  the  responsibility  of  designers  to  ensure 
that  electronic  components  operate  efficiently  and  effectively 
throughout  the  specified  thermal  limits.  Therefore,  it  is 
extremely  important  that  design  engineers  have  the  capability 
to  accurately  and  rapidly  predict  the  temperature  distribution 
on  multilayered  structures  prior  to  prototype  production.  The 
overriding  reasons  for  performing  a  precise  thermal  analysis 
are  to  increase  component  reliability,  ensure  proper  material 
selection,  ensure  bias  stabilization,  and  reduce  or  eliminate 
the  possibility  of  catas'*'rophic  thermal  failure .  [Ref .  1] 
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A.  RELIABILITY 

There  is  a  predictable  relationship  between  the  operating 
temperature  of  electronic  components  and  reliability  [Ref.  2]. 
The  materials  used  in  the  fabrication  of  components  have 
temperature  limitations.  Should  these  temperature  boundaries 
be  exceeded,  the  physical  and  chemical  properties  of  the 
material  are  altered  and  the  device  fails.  Figure  3  displays 
the  intimate  relationship  between  failure  rate  and  component 
operating  temperature  for  some  selected  devices.  Furthermore, 
it  is  an  established  fact  that  the  reliability  of  an 
electronic  component  is  inversely  proportional  to  the  junction 
or  component  temperature  and  is  also  directly  linked  to 
failure  rates  [Ref.  1]. 

Consider  Figure  4  which  illustrates  the  "bathtub" 
mortality  curve  with  the  failure  rate  of  a  particular 
component  plotted  against  component  age  during  operation 
within  thermal  limits.  The  high  failure  rate  in  the  interval 
prior  to  t^,  also  known  as  the  burn-in  period,  is  considered 
to  be  the  result  of  poor  quality  control  during  the 
fabrication  process .[ Ref .  1] 

The  area  of  highest  concern  is  the  interval  between  t^  and 
t^.  This  period  is  considered  to  be  the  useful  life,  since 
with  proper  quality  control,  testing  and  burn-in  procedures, 
t^  is  equal  to  zero.  Failures  that  occur  in  this  interval  are 
due  to  a  variety  of  causes  and  are  unpredictable .[ Ref .  1] 
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Figure  3.  Failure  rate  vs.  Temperature  for 
selecte  devices  [Ref.  2]. 
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Failure  also  occurs  gradually  due  to  sustained  operations 
within  specified  temperature  limits.  If  all  aspects  of 
fabrication  have  been  performed  correctly  only  a  small 
fraction  of  components  will  have  failed  prior  to  reaching  t^, 
also  known  as  the  wear-out  period.  Failures  that  occur  in  this 
period  are  due  to  the  slow  and  never  ending  deterioration  of 
materials . 

B.  MATERIAL  SELECTION 

The  fabrication  of  electronic  components  results  in  the 
joining  of  several  different  materials.  Consider  Figure  2 
which  depicts  a  typical  semiconductor  structure.  When  power  is 
applied  to  the  components,  heat  is  dissipated  to  the  substrate 
and  subsequently  to  the  carrier.  The  mechanical  properties  of 
these  materials  are  all  affected  differently  by  changes  in 
temperature.  Opposing  mechanical  and  chemical  reactions  due  to 
environmental  conditions  and  contaminants  may  result  in 
component  performance  degradation  or  a  reduction  in  useful 
life.  Table  1  lists  temperature  related  factors  that  may 
affect  component  performance. 

The  primary  objective  in  the  selection  of  materials  for 
the  fabrication  of  an  electronic  assembly  is  to  achieve  the 
desired  level  of  correlation  within  the  finished  product.  As 
packaging  densities  increase  thermal,  mechanical,  electrical, 
and  chemical  coupling  becomes  very  strong.  This  high  level  of 
coupling  can  be  both  an  advantage  and  a  disadvantage.  For 
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example,  a  high  level  of  correlation  is  desirable  during 
fabrication  to  ensure  an  uniform  product.  However,  in  use, 
strong  coupling  is  generally  more  desirable  for  moderate 
temperature  deviations  and  weaker  coupling  is  more  desirable 
for  large  temperature  deviations.  In  the  case  of  large 
temperature  deviations  strong  coupling  may  result  in  the 
catastrophic  failure  of  many  connected  components  while  weak 
coupling  may  limit  the  number  of  failed  components.  Therefore, 
the  strength  of  coupling  between  materials  must  be  based  on 
the  type  of  failure  most  likely  to  occur  and  an  accurate 
thermal  analysis  must  supply  this  information  to  assist  in 
proper  material  selection  [Ref.  3]. 


TABLE  1.  TEMPERATURE  FACTORS  [Ref.  4] 


Mechanism 

Effect  on  Equipment 

Accelerating  factors 

Increasing 

T  emperature 

Loss  of  strength,  reduced 
stiffness,  reduced  resonant 
frequency,  softening,  distortion, 
aging,  and  creep 

Lubricants,  rubber  parts, 
plastics,  corrosion,  fatigue, 
load  intensity,  and  time 
durat i on 

Reducing  femperature 

Increased  viscosity,  increased 
stiffness,  increased  resonant 
frequency,  brittleness,  and 
reduced  impact  resistance 

Lubricants,  rubber  parts, 
plastics,  and  time  duration 

Thermal  Expansion  and 
Contraction 

Change  in  size  and  shape, 
buckling,  cracking,  distortion, 
and  loosening 

Temperature  cycling,  temperature 
range,  unequal  expansion 
coefficients,  stress 
concentrations 
,  and  lack  of  strain  relief 
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C.  BIAS  STABILIZATION 

The  first  step  in  the  design  and  implementation  of  a 
semiconductor  device  is  to  establish  a  stable  and  predictable 
electrical  operating  point.  This  procedure,  known  as  bias 
stabilization,  attempts  to  determine  a  stable  operating  point 
that  is  virtually  independent  of  external  component 
parameters.  However,  as  external  parameters  change,  the 
operating  point  is  directly  affected.  Therefore,  a  good  bias 
design  ensures  that  components  will  always  operate  within  a 
certain  range  of  their  nominal  value. [Ref.  5] 

Consider  Figure  5  which  displays  a  transistor  connected  in 
the  common-emitter  configuration.  Suppose  that  a  proposed 
operation  requires  a  specific  collector  to  emitter  voltage 
(Vj,g)  .  The  circuit  consists  of  a  battery  or  some  other  source 
that  provides  a  bias  voltage  the  collector  resistor  R^, 
and  the  transistor.  By  Kirchoff's  voltage  law 

-V^cc  ^  *  ^CE  =  0 


which  results  in  a  collector  to  emitter  voltage  of 

Should  the  collector  current  be  allowed  to  increase  in  excess 
of  tolerable  limits,  must  decrease  because  and  R^.  are 
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Figure  5.  Schematic  of  transistor  connected  in  common-emitter 
configuration  [Ref.  1]. 
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fixed  values.  Therefore,  it  should  be  noted  that,  if  a  high 


junction  temperature  causes  an  increase  in  1^,,  can  no 
longer  be  maintained  at  the  level  desired  to  perform  the 
desired  operation. [Ref .  1] 

1.  Operating  in  the  Forward  Bias  Region 

As  an  example  of  electronic  component  temperature 
dependence,  consider  a  diode  operating  in  the  forward  bias 
region.  In  the  forward  region  the  i-v  relationship  is  closely 
approximated  by 


V 

i  =  15(0"^-  -  1) 


(3) 


In  this  equation  I5  is  a  constant  for  a  given  diode  at  a  given 
temperature.  The  current  Ij  is  usually  called  the  saturation 
current.  However,  another  name  for  it  is  the  scale  current, 
which  arises  from  the  fact  that  I^  is  directly  proportional  to 
the  cross-sectional  area  of  the  diode.  Furthermore,  it  can  be 
seen  in  Table  2  that,  I5  is  a  very  strong  function  of 
temperature , [Ref .  5] 

The  temperature  relationship  between  l^  and  the 
forward  current  i  is  derived  from  the  voltage  V^.  This 
constant,  called  the  thermal  voltage,  is  given  by 

Vj.  =  —  (4) 

Q 
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where 


k  =  Boltzman's  constant,  1.38- 10'^^  J/  K 
T  =  the  absolute  temperature,  K 
q  =  the  charge  on  the  electron,  1.602- 10'^’  C 

Table  2  illustrates  this  relationship  and  emphasizes  the  need 
to  accurately  analyze  a  proposed  assembly  prior  to 
fabrication. 


TABLE  2.  TEMPERATURE  DEPENDENCY  OF  Ig  ON  i 
FOR  SELECTED  MATERIALS 


GERMANIUM 

SILICON 

“C 

i 

Is 

i 

25 

3.0  /iA 

18.01  juA 

50.0  ryA 

82.31  rjh 

95 

0.384  mA 

1.473  mA 

51.2  mA 

61.39  mA 

165 

49 . 2  itiA 

0.136  A 

52.4  mA 

49.2  mA 

D.  CATASTROPHIC  THERMAL  FAILURE 

Another  of  the  primary  goals  of  techniques  in  advanced 
thermal  control  is  to  provide  a  thermal  environment  for  a 
diversity  of  components  that  are  in  increasingly  close 
proximity  to  each  other.  Figure  6  illustrates  the  increasing 
level  of  packaging  densities.  With  increasing  complexity  comes 
an  increased  level  of  connections  between  dissimilar  material 
and  a  greater  possibility  for  exceeding  temperature 
limitations.  Therefore,  it  is  necessary  that  designers  have 
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Components  per  circuit 
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some  knowledge  of  the  possible  range  of  component  and 
environmental  variations  in  order  to  prevent  catastrophic 
thermal  failure . [Ref .  1] 

Catastrophic  thermal  failure  is  defined  as  an  immediate, 
thermally  induced,  total  loss  of  electronic  function  in  a 
specified  component.  This  type  of  failure  is  the  result  of  a 
component  melting  due  to  excessive  temperature,  a  thermal 
fracture  of  the  substrate  or  carrier,  or  a  separation  of  leads 
and  the  external  network.  It  is  generally  considered  to  be 
dependent  on  the  local  temperature  field,  operating  history, 
and  operating  modes  of  the  component.  As  previously  stated,  a 
variety  of  problems  arise  when  components  are  subjected  to 
temperatures  in  excess  of  their  rated  limits.  Furthermore,  it 
is  extremely  difficult  to  determine  the  precise  temperature  at 
which  catastrophic  failure  may  occur.  The  incorporation  of  an 
accurate  thermal  analysis,  in  combination  with  test  and 
operating  experience,  may  be  used  to  generate  a  catastrophe 
free  upper  operating  limit.  These  maximum  allowable  operating 
temperatures  are  used  to  generate  the  master  thermal  control 
configuration  for  the  system . [Ref . 1 ] 
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III. 


HEAT  TRANSFER 


Heat  transfer  is  defined  as  all  energy  flows  that  arise  as 
a  result  of  temperature  differences  [Ref.  6].  Because 
electronic  components  are  not  one  hundred  percent  efficient, 
they  produce  heat  as  well  as  the  desired  output.  In  the  case 
of  semiconductor  devices,  heat  develops  in  parts  having  low 
thermal  efficiencies,  such  as  the  die.  One  of  the  major 
objectives  of  packaging  is  to  develop  an  effective  system  for 
the  removal  of  heat  from  these  parts  [Ref.  4].  It  is 
imperative  that  design  engineers  understand  all  modes  of  heat 
transfer  in  order  to  incorporate  an  efficient  method  of  heat 
removal  into  component  designs.  The  modes  of  heat  transfer  are 
conduction,  convection,  and  radiation. 

A.  CONDUCTION 

Conduction  is  the  transfer  by  molecular  motion  of  heat 
between  one  part  of  a  body  to  another  part  of  the  same  body  or 
between  one  body  and  another  in  physical  contact  [Ref.  1]. 
Joseph  Fourier,  a  French  physicist,  proposed  that  the  rate  of 
heat  flow  through  a  material  by  conduction  is  proportional  to 
the  area  o:^  the  material  normal  to  the  heat  flow  path  and  to 
the  temperature  gradient  along  the  heat  flow  path. 
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This  proportionality  is  represented  mathematically  by 


q  « 


(5) 


where  the  minus  sign  allows  for  a  positive  heat  flow  in  the 
presence  of  a  negative  temperature  gradient.  The  introduction 
of  a  proportionality  constant,  known  as  thermal  conductivity, 
results  in  the  following  rate  equation  which  describes  this 
mechanism  [Ref,  1]: 


Q  = 


dx 


(6) 


where 

k  =  thermal  conductivity  of  the  material,  W/m-  C 
A  =  area  of  the  heat  flow  path,  m^ 

dT/dx  =  change  in  temperature  per  unit  length.  C/m 
q  =  rate  of  heat  flow,  W 

1.  General  Equation  of  Heat  Conduction 

The  first  step  in  the  analytical  solution  of  a  heat 
conduction  problem  for  a  given  structure  is  to  choose  an 
orthogonal  coordinate  system  such  that  the  surfaces  coincide 
with  the  boundary  surfaces  of  the  structure  [Ref.  7].  In  the 
case  of  the  model  builder  developed  in  this  thesis, the 
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rectangular  coordinate  system  will  be  employed.  The  general 
equation  of  heat  conduction  is  given  as 


Then,  assuming  k,  C,  and  p  are  independent  of  temperature, 
direction,  and  time,  the  resulting  equation  is 

d^T  ^  ^  ^_T  ^  q  ^  2^^  .q. 

dx‘^  dy^  dz^  ^  a 


where 


T  =  temperature,  C 

X,  y,  and  z  =  cartesian  coordinates,  m 
t  =  time,  sec 

k  =  thermal  conductivity,  W/m-  C 
q  =  internal  heat  generation,  W/m^ 
a  =  thermal  diffusivity,  k/pC,  mVsec 

There  are  several  variations  of  the  general  equation 
of  conduction.  The  first,  known  as  the  Fourier  equation, 
provides  a  solution  for  a  system  that  contains  no  heat 
sources : 

d^T  ^  d^T  (9) 

dx^  dy^  dz^  a  8t 
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The  second  variation,  known  as  the  Poisson  equation, 
supplies  a  solution  for  a  system  in  which  the  temperature  does 
not  vary  with  time: 


dx^  dy^  dz^  k 


(10) 


The  third  and  final  variation  of  the  general  equation 
of  conduction  provides  a  solution  for  a  system  void  of  heat 
sources  and  operating  in  steady  state.  The  resulting  equation, 
known  as  the  Laplace  equation,  is  given  as 


^  .  Zl  =  v=r  =  0 

dx^  dy^  dz^ 


(11) 


2.  Simple  Plane  Slab 

Consider  Figure  7  which  illustrates  a  simple  plane 
slab  with  face  temperatures  T,  and  T^ .  Using  only  one 
dimension,  equation  (11)  is  reduced  to 

T 

-T-J.  =  0 
dx^ 


By  integrating  twice  and  applying  boundary  conditions  the 
temperature  distribution  across  the  slab  is  seen  to  be 

T  =  T.  -  —  (T,  -  T.)  (12) 

.  L 
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Figure  7.  Conduction  heat  transfer  through  a  simple 
plane  slab. 
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Insertion  of  equation  (12)  into  equation  (6)  produces  a 
solution  for  the  heat  flow  across  the  slab: 


g  = 


(13) 


Ohm's  Law  indicates  a  direct  analoqy  between  heat 
flow,  equation  (13),  and  the  flow  of  electrical  current 
through  a  resistor,  V  =  RI .  This  electrothermal  analogy  is 
extremely  useful  in  the  solution  of  one  dimensional,  steady 
state  problems  without  energy  generation  and  will  be  developed 
further  in  the  next  section. [Ref .  8] 

3.  Electrothermal  Analog 

As  previously  stated,  there  is  a  direct  analogy 
between  heat  flow  across  a  simple  plane  slab,  equation  (13) , 
and  electrical  current  governed  by  Ohm's  law: 


I  = 


V 

R 


(14) 


In  this  case,  the  analogous  quantities  are 

Current  I  «  Heat  Flow  q 
Potential  V  «  Temperature  Difference  AT 
Resistance  R  Thermal  Resistance  R 
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It  is  easily  seen  that  for  the  heat  flow  in  a  simple 
plane  slab  described  by,  equation  (13),  the  thermal  resistance 

is 


R  = 


AT 

Q 


L 

kA 


(15) 


The  electrothermal  analog  for  conduction  across  a  simple  plane 
slab  is  shown  in  Figure  8. [Ref.  1] 


B .  CONVECTION 

Convection  is  defined  as  the  process  by  which  thermal 
energy  is  transferred  to  or  from  a  solid  by  a  fluid  flowing 
past  it.  Should  the  fluid  flow  be  the  result  of  a  temperature 
difference  the  phenomena  is  called  natural  or  free  convection. 
On  the  other  hand,  when  a  pump  or  fan  causes  the  mass  movement 
the  process  is  called  forced  convection. [Ref.  1] 

Recall  that  at  the  interface  between  a  solid  and  a  fluiJ 
that  heat  is  transferred  by  conduction  and  must  obey  Fourier's 
law,  equation  (6).  Due  to  the  difficulty  encountered  in 
accurately  measuring  the  temperature  gradient,  Newton 
suggested  that  the  surface  heat  transfer  rate  be  related  to 
the  product  of  surface  area  and  the  temperature  difference 
between  the  surface  and  the  fluid.  The  results  of  this 
proposition  lead  to  Newton's  law  of  cooling: 


q  =  hA{T^  -  Tf) 


(16) 
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where  h  is  a  proportionality  factor  that  has  become  known  as 
the  surface  heat  transfer  coefficient. [Ref .  1] 

From  a  comparison  of  Newton's  law  of  cooling,  equation 
(16),  and  Fourier's  law,  equation  (6),  one  can  derive  that  the 
surface  heat  transfer  coefficient  can  be  related  to  the 
thermal  conductivity,  the  wall  temperature  gradient  of  the 
fluid,  and  the  surface  fluid  temperature  difference: 

h  -  -3^  = 

A■^T  AT 


Therefore,  any  correlation  between  heat  transfer  coefficients 
must  reflect  the  dependence  of  h  on  the  thermal  conductivity 
of  the  fluid  and  on  the  ratio  of  the  wall  temperature  gradient 
to  the  temperature  difference. [Ref .  1] 

1.  Electrothermal  Analog 

The  addition  of  heat  transfer  by  convection  to  both 
surfaces  of  the  simple  plane  slab  of  Figure  7  results  in  the 
configuration  shown  in  Figure  9.  In  the  convective  case 
thermal  resistance  is  represented  by 


R  = 


hA 


(18) 
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^  _ _ _  *  -  .  ^ 


Figure  9.  Convective  heat  transfer  on  a  simple  plane  slab. 
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Therefore,  the  total  thermal  resistance  is 


R  = 


1 

h^A 


L  ^  1  ^  2  J_ 

kA  h^A  A 


(19) 


which  is  represented  by  the  electrothermal  analog  shown  in 
Figure  10.  A  simple  consideration  of  circuit  theory  then  shows 
that: 


AT 

(T,  -  Tj) 

R 

1 

1 

A 

h,  k 

^2. 

C.  RADIATION 

Heat  transfer  by  radiation  is  the  means  by  which  thermal 
energy  can  be  transmitted  through  a  space  without  an 
intervening  medium  while  obeying  the  laws  of  electromagnetics. 
Thermal  radiation,  while  traveling  at  the  speed  of  light,  may 
be  absorbed,  reflected,  or  transmitted  upon  contact  with  a 
surface.  An  ideal  black  body  absorbs  all  incident  radiation 
and  reflects  and  transmits  none  of  it.  The  concept  of  the 
black  body  is  useful  because  laws  governing  its  radiation  are 
simple  and  many  real  bodies  may  be  treated  approximately  as 
black  bodies  [Ref.  1]. 

Materials  used  in  the  fabrication  of  electronic  components 
are  classified  as  gray.  Gray  bodies  are  diffusely  reflecting 
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V  V  V 


Figure  10.  Conduction  and  convection  electrothermal 
equivalent . 
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opaque  surfaces  [Ref.  4].  These  surfaces  reflect  equal  amounts 
of  energy  over  the  thermal  radiation  spectrum  (wavelengths  of 
about  0.1  /im  to  about  100  ixm)  in  all  directions.  The  heat 
transfer  efficiency  of  this  mode  depends  on  the  configuration, 
the  orientation,  and  the  temperatures  of  the  surfaces  in  the 
electronic  assembly. 

1.  Transformation  of  the  General  Radiation  Equation 

The  use  of  the  thermal  radiation  equation  in 
analytical  studies  is  made  difficult  by  its  dependence  on  the 
fourth  power  relationship  between  the  temperatures.  Due  to  the 
nonlinear  characteristics  and  the  complexity  of  the 
calculations  a  computer  program  is  the  desired  method  to  solve 
problems  that  have  a  significant  transfer  by  radiation.  The 
general  equation  for  radiation  interchange  is: 

q  =  oF^F^  (Tt  -  Tt)  (21) 

where 

a  =  Stefan-Boltzman  constant,  5.669- 10  ®  W/m^- 

=  shape  factor  that  accounts  for  the 
arrangement  of  the  of  the  radiating  source  and 
absorbing  receiver 

=  emissivity  factor  that  accounts  for  the 
ability  of  the  source  and  receiver  to  emit  or  absorb 
radiation 

T^  =  temperature  of  the  source,  K 
T^  =  temperature  of  the  receiver,  K 
A  =  surface  area,  m^ 
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It  is  important  to  note  that  the  absolute  temperature  scale 
must  be  used  when  considering  radiation. 

One  way  to  handle  computations  involving  heat  transfer 
by  radiation  is  to  transform  the  general  radiation  equation, 
equation  (21),  into  a  configuration  compatible  with  Fourier's 
law.  In  this  case,  linearization  of  the  general  radiation 
equation  is  the  method  used  to  produce  the  desired  result. 
This  is  achieved  by  factoring  the  difference  in  the  fourth 
power  of  the  temperatures  as  follows: 

[ji  -  Ti)  =  VIZ  -  T~)  (I'i  -  r;) 

=  {T;  *  li)  (T^  T^)  (Tg  -  T^)  (22) 

Inserting  this  into  equation  (21)  results  in 

g  =  aF^F^A  ( +  2^)  ( (23) 


A  radiative  heat  transfer  coefficient  may  therefore  be  defined 
as 

h,  =  oF^,F^(TI  *  li)  (T^  ^  T,)  (24a) 

or 

=  oF^F^ili  -  TiT^  ^  T^tI  +  Ti)  (24b) 
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Then,  substituting  h^.  into  equation  (23)  ,  radiation  heat 
transfer  may  be  treated  exactly  as  convection  at  the  boundary. 

A  thermal  resistance  for  radiation  heat  transfer  can 
now  be  proposed: 


Figure  11  provides  an  illustration  of  the  electrothermal 
equivalent  with  the  addition  of  radiation  resistance  in 
parallel  with  convective  resistance. 


Figure  11.  Radiation,  convection,  and  conduction 
electrothermal  equivalent. 
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2.  General  Problems  of  Heat  Transfer  by  Radiation 

In  the  calculation  of  heat  transfer  by  radiation,  it 
is  usually  necessary  to  approximate  real-body  behavior  by  the 
gray-body  idealization.  The  assumption  that  emissivity  is 
always  equal  to  absorptivity  is  frequently  required  if  the 
problem  is  to  be  solved.  Even  a  simple  situation  becomes  quite 
complex  if  real  body  behavior  is  considered  and  the  resulting 
convenience  of  replacing  absorptivity  with  emissivity  is  lost. 
A  second  difficulty  in  considering  real-body  behavior  is  the 
lack  of  sufficient  data.  To  properly  account  for  real-body 
behavior,  extensive  tabulations  would  be  required. [Ref .  6] 

The  usual  problem  of  heat  transfer  by  radiation  is 
further  complicated  not  only  because  the  surfaces  are  nonblack 
but  also  because  the  configuration  of  the  areas  involved  is 
not  simple  and  reflecting  surfaces  may  be  present  to  augment 
the  direct  exchange  [Ref.  6].  There  are  a  number  of  methods 
available  to  assist  in  the  solution  of  heat  transfer  by 
radiation,  however,  their  development  is  beyond  the  scope  of 
this  thesis. 
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IV. 


FINITE  DIFFERENCES 


The  need  to  utilize  a  computer  to  determine  the 
temperature  distribution  within  an  electronic  component  or 
system  has  made  finite  difference  methods  very  desirable.  This 
numerical  method  is  useful  in  problems  involving 
nonlinearities,  complex  geometries,  complicated  boundary 
conditions,  or  a  system  of  coupled  partial  differential 
equations.  The  purpose  of  this  section  is  to  provide  the 
reader  with  some  basic  concepts  involved  in  finite  difference 
methods  for  solving  differential  equations.  Furthermore,  it 
demonstrates  the  methodology  used  to  formulate  n-node 
equations  in  n-unknown  temperatures  and  instills  confidence 
that  this  numerical  method  is  capable  of  generating  an 
accurate  thermal  model. 

A .  FUNDAMENTAL  CONCEPTS 

Consider  equation  (7),  the  general  equation  for  heat 
transfer  by  conduction.  In  order  to  produce  a  numerical 
solution  to  a  conduction  heat  transfer  problem  it  is  necessary 
to  reconfigure  the  partial  differential  equation  into  a  form 
that  allows  differentiation  to  be  performed  by  numerical 
methods.  Therefore,  it  is  essential  that  accurate 
approximations  of  the  first  and  second  derivatives  be 
obtained . 
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1.  First  Derivative  Approximation 

The  derivative  of  a  function  at  any  point  can  be 
expressed  by  a  finite  difference  approximation  by 
incorporating  a  Taylor  series  expansion  about  that  point. 
Consider  Figure  12  where  T(x)  is  a  function  that  can  be 
expanded  by  a  Taylor  series.  The  Taylor  series  expansions  of 
the  functions  T(x+Ax)  and  T(x-Ax)  about  a  point  .v  are: 

T(x*£i,x)  =  r(x)  »  Axr(x)  +  A^^^t"(x)  *  iA^rW(x)  (26a) 

r(x-Ax)  =  7(x)  -  AxT^x)  +  1^^t"(x)  -  JA^r'"(x)  (26b) 

In  order  to  determine  the  first  derivative,  equations  (26a) 
and  (26b)  are  solved  for  T'(x). 

{ x)  =  r(x  Ax)  ~  7(x)  _  (Ax')  ^/f  ^  __  ( Ax)  j.///  ^  _  (27) 


'(X) 


Tjx)  -  T(x  -  Ax; 

Ax 


(Ax)  _  (Ax) 

2  6 


r'"(x) 


(28) 


Subtracting  equation  (26b)  from  (26a)  and  solving  for  T' (x) 
produces 


T'(x) 


T(x  Ax)  -  T(x  -  Ax) 
2Ax 


( Ax)  ^ 
6 


T'"  {x) 


(29) 
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From  equations  (27),  (28),  and  (29)  the  first 

derivative  approximations  of  T' (x)  are  determined  to  be 


T'{x)  =■  —  =  "  Tjx) 


forward  difference 


T'ix)  =  ^  = 


dT  _  T{x)  -  T(x  -  Ax) 


backward  difference 


T'ix)  = 


dT  _  T(x  *  Ax)  -  T(x  -  Ax) 


central  difference 


and  for  all  it  is  observed  that 


^  ^ 

dx  Ax 

2.  Second  Derivative  Approximation 

To  obtain  the  second  derivative  of  the  function  T(x) , 
Figure  12,  consider  the  Taylor  series  expansions  of  the 
functions  T(x+2Ax)  and  T(x-2Ax)  about  a  point  x. 


T{x*2Ax)  =  Tix)  ^  lAxT'ix)  *  2(.Ax)^T"(x)  +  (AX)  ^r"'(x)  +...  (33a) 


T(x-2Ax)  =  Tix)  -  2AxT'ix)  *  2iAx)^T"{x)  -  ( Ax) ’T'" (x) 


(33b) 
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Inserting  the  corresponding  first  derivative  approximation 
into  equations  (33a)  and  (33b)  and  solving  for  T" (x)  results 
in 

r'(x)  ^  ^  ^  ^  ~  foiward  difference  (34) 

dx^  (Ax)^ 


-"(x)  =  ^  =  r(x  -  2Ax)  ^  T(x)  -  2Tix  :  Ax)  tackward  difference  (35) 
dx^  iAx)^ 

The  central  difference  is  obtained  by  eliminating  T'  (x) 
between  equations  (26a)  and  (26b) . 

T‘'(x)  =  -2-T  =  ~  *  T(x  >  Ax)  -  2T(x)  central  difference  (36) 

dx^  (Ax)^ 

and  for  all  it  is  observed  that 

T»(x)  .  ^ 

dx^  (Ax)^ 


B.  NODE  ANALYSIS 

The  first  step  in  the  physical  formulation  of  a  solution 
by  node  analysis  is,  as  previously  stated,  to  divide  the 
region  into  a  finite  number  of  subvolumes  as  shown  in 
Figure  13.  The  centroid  of  each  subvolume  is  called  a  node  and 
is  considered  to  be  representative  of  the  entire  subvolume. 
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Figure  13.  Structure  divided  into 
subvolumes . 
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At  this  point,  continuity  must  be  applied.  In  the  steady 
state,  the  algebraic  sum  of  all  the  heat  leaving  any  node  must 
equal  zero.  As  an  example,  consider  an  interior  node  on  the 
upper  layer;  select  node  number  14.  Figure  14  provides  an 
expanded  picture  around  node-14.  In  order  to  maintain  clarity 
the  external  environment  will  be  considered  as  a  single  node, 
number  81.  The  resulting  energy  balance  or  node  equation  can 
be  written  as 

^14.10  ^  ^14,1^  ^  ^14,15  ^  ^14.33  ^  ^14 , 31  +  =  0  (37) 


Here,  each  numerical  subscript  represents  the  heat  flow  from 
node-14  to  the  indicated  node.  Furthermore,  allows  for 
energy  input  by  an  external  source  such  as  a  dissipating 
component . 

The  energy  balance  simply  supplies  the  pertinent  energy 
terms.  It  is  important  to  note  that  differences  in  thicknesses 
and  thermal  conductivities  must  be  taken  into  consideration. 
The  rate  equations  may  be  written  in  terms  of  thermal 
resistances  and  temperature  differences: 


ic.AxAz, 

^14.13  =  -(^4  -  ^13) 


(38a) 


k.AxAz. 

<714,15  =  —  2^^  -(^14  -  T’xs) 


(38b) 
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(38c) 


^14, 10 


iCj^AyA 

Ax 


(7’,,  -  T,,) 


Q. 


4 . : 


k^Ay^z^ 

Ax 


(T., 


T-.e) 


(38d) 


'^14 ,61 


2k^AxAy 

'^1 


(^14  - 


‘  81  ‘ 


(38e) 


»14 .38 


2  AxAy 


(4^) 


( 


Azo 


( T 

'  14 


^38) 


(38f) 


where 


and  k,  =  thermal  conductivities  of  the  top  layer 
and  middle  layer,  respectively,  W/m-  C 

Ax  and  Ay  =  path  lengths  in  their  respective 
direction  and  are  node  dimensions  for  surface 
area  calculations,  m 

Az,  and  AZj  =  the  thicknesses  for  the  upper  and 
middle  layers,  respectively,  m 

=  temperature  of  the  node  indicated  by  the 
subscript,  C 


Substituting  these  rate  equations  into  the  energy  balance 
equation,  equation  (37),  results  in  an  equation  containing 
five  unknown  temperatures. 
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For  simplicity  let  Ax=Ay=Az^=A22=l  and  k^=k2=k.  Then 


kcr,,  -  T,,)  -  k(T,,  -  r,3)  -  k(T,,  -  T-s)  -  k(T,,  -  T, 


13' 


2k(r,,  -  Tg.)  +  k(r,,  -  Tjs)  +  g,  =  0 


(39) 


further  reduction  results  in 


-T. 


10 


•  13 


T  -  T 
-‘is  -‘la 


7r,,  -  T. 


‘  14 
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Si  ^  2  T, 


81 


(40) 


An  equation  of  this  type  can  be  written  for  every  node  in  the 
region  whose  temperature  is  unknown.  This  results  in  a  set  of 
n-equations  that  relate  the  /j-unknown  temperatures.  There  are 
a  variety  of  methods  available  to  solve  for  the  unknown 
temperatures.  In  large  structures  the  most  desirable  method 
would  be  a  thermal  analysis  program. 
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V. 


THE  MODEL  BUILDER 


As  previously  stated,  the  model  builder  currently  being 
used  by  the  thermal  analysis  software  package  requires  an 
inordinate  amount  of  time  for  data  entry.  This  thesis  presents 
a  specific  model  builder,  TMDL,  that  will  become  part  of  the 
software  package.  TMDL  provides  a  user  friendly,  rapid  method 
to  assist  in  the  development  of  a  thermal  model  with  a 
specific  geometry  as  in  Figure  2.  Specifically,  it  generates 
a  properly  formatted  output  data  file  for  use  by  the  thermal 
analyzer  enroute  to  preparing  an  accurate  listing  of  the 
temperature  distribution  of  the  structure. 

A.  THE  THERMAL  ANALYZER  INPUT  DATA  FILE 

The  model  builder  generates  a  data  file  from  physical 
characteristics  of  the  structure  provided  by  the  user.  This 
data  file,  also  called  the  thermal  analyzer  input  data  file, 
must  be  in  a  format  that  is  completely  acceptable  to  the 
thermal  analyzer.  Furthermore,  it  will  consist  of  five  lines 
and  up  to  seven  data  sets.  This  section  describes  each  line 
and  data  set  and  their  relationship  with  TMDL.  Figure  15 
displays  a  partial  data  file. 

Line  one  is  the  title  line.  It  may  be  left  blank  or  may 
contain  up  to  79  alphanumeric  characters.  The  user  selected 
title  appears  at  the  top  of  the  data  file. 


EXAMPLE  OUTPUT 

DATA  FILE 

740 

6 

0 

0  0 

0 

0  0 

2 

0 

0 

0 

750  50 

A 

6 

2  4 

6 

0  0 

0 

.0500000 

. 6666700 

12 

.8000000  78. 

00000 

76.500000000 

76.900000000  75.900000000 

77.100000000 

1  77.400000000  75.6000 

7  7551 

21 

7521 

111 

7511 

2411 

73.700 

36.850 

87 .709 

43.855 

2 . 047 

1 .456 

6  11 

31 

7521 

121 

7511 

2421 

36.850 

36.850 

87.709 

43.855 

2.047 

1.456 

7  21 

41 

7521 

131 

7511 

2431 

36.850 

36.850 

87.709 

43.855 

2.047 

1.456 

6  31 

51 

7521 

141 

7511 

2441 

36.850 

36.850 

87.709 

43.855 

2.047 

1.456 

6  41 

61 

7521 

151 

7511 

2451 

36.850 

36.850 

87.709 

43.855 

2.047 

1.456 

6  51 

71 

7521 

161 

7511 

2461 

36.850 

36.850 

87.709 

43.855 

2.047 

1.456 

7  61 

81 

7521 

171 

7511 

2471 

36.850 

36.850 

87.709 

43.855 

2.047 

1.456 

6  71 

91 

7521 

181 

7511 

2481 

36.850 

36.850 

87.709 

43.855 

2.047 

1.456 

6  81 

101 

7521 

191 

7511 

2491 

36.850 

36.850 

87.709 

43.855 

2.047 

1.456 

7  91 

7541 

7521 

201 

7511 

2501 

36.850 

73.700 

87.709 

43.855 

2.047 

1.456 

6  7551 

121 

11 

211 

7511 

2511 

73.700 

36.850 

43.855 

43.855 

2.047 

1.456 

6  111 

131 

21 

221 

7511 

2521 

36.850 

36.850 

43.855 

43.855 

2.047 

1.456 

7  121 

141 

31 

231 

7511 

2531 

36.850 

36.850 

43.855 

43.855 

2.047 

1.456 

6  131 

151 

41 

241 

7511 

2541 

36.850 

36.850 

43.855 

43.855 

2.047 

1.456 

6  141 

161 

51 

251 

7511 

2551 

36.850 

36.850 

43.855 

43.855 

2.047 

1.456 

6  151 

171 

61 

261 

7511 

2561 

36.850 

36.850 

43.855 

43.855 

2.047 

1.456 

6  161 

181 

71 

271 

7511 

2571 

36.850 

36.850 

43.855 

43.855 

2.047 

1.456 

6  171 

191 

81 

281 

7511 

2581 

36.850 

36.850 

43.855 

43.855 

2.047 

1.456 

7  181 

201 

91 

291 

7511 

2591 

36.850 

36.850 

43 . 855 

43.855 

2 . 047 

1.456 

6  191 

7541 

101 

301 

7511 

2601 

36.850 

73.700 

43.855 

43.855 

2.047 

1.456 

6  7551 

221 

111 

311 

7511 

2611 

73.700 

36.850 

43.855 

43.855 

2 . 047 

1.456 

6  211 

231 

121 

321 

7511 

2621 

36.850 

36.850 

43.855 

43.855 

2 . 047 

1.456 

Figure  15.  Output  data  file. 
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Line  two  is  the  problem  data  line.  It  has  nine  entries  of 
which  two  are  under  user  control,  number  of  nodes  under 
consideration  and  unit  type.  One  entry,  the  number  of  constant 
temperatures  is  preset  at  six  for  this  specific  model.  The 
remaining  entries  have  applications  to  models  associated  with 
heaters,  unique  exponents,  secondary  heat  input,  temperature 
coefficients  and  curves,  and  nodes  controlling  fast  heat. 
These  entries  are  not  applicable  to  this  model  and  are  preset 
to  zero. 

Line  three  places  a  zero  at  three  points  and  is  beyond  the 
user's  control.  Therefore,  no  further  discussion  is  required. 

Line  four  is  the  problem  capability  line.  This  line 
defines  the  maximum  values  for  the  entries  in  line  two.  The 
first  entry  is  750,  the  number  of  nodes  for  which  the  analysis 
is  dimensioned.  This  is  significant  because  the  first  constant 
temperature  will  be  assigned  the  node  number  751.  The  second 
entry  is  50  which  is  the  maximum  number  of  constant 
temperatures  in  accordance  with  the  analyzer  dimension 
statement.  Therefore,  it  is  possible  to  have  50  constant 
temperatures  allocated  from  node  751  to  800.  The  third  entry 
is  preset  to  six  for  an  application  concerning  heaters  and  is 
not  applicable  to  this  model.  The  balance  of  the  entries  in 
line  four  represent  a  listing  of  data  sets  that  are  required 
for  the  particular  analysis  at  hand.  TMDL  uses  three  data  sets 
that  will  be  discussed  in  what  follows. 
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Line  five  contains  five  items  that  concern  the  level  of 


accuracy  that  the  thermal  analyzer  will  achieve.  These  entries 
are  preset.  Item  one  provides  the  desired  level  of  accuracy 
between  iterations.  Because  the  thermal  analyzer  solves 
iteratively,  there  must  be  an  error  criterion  at  which  point 
calculations  cease.  A  number  that  is  too  small  will  cause  the 
computer  to  run  excessively  and  a  number  that  is  too  large 
will  not  give  the  desired  accuracy.  Therefore,  a  trade  off  is 
necessary  and  such  a  tradeoff  has  indicated  that  C.05  is 
satisfactory.  Item  two  is  a  damping  factor  that  is  used 
between  iterations  to  prevent  temperature  oscillations  between 
iterations.  The  third  entry  is  the  maximum  number  of 
iterations.  This  prevents  excessive  computer  time  in  the  event 
of  faulty  input  data.  Item  four  is  a  convergence  factor  that 
adjusts  the  damping  to  close  to  the  critical  value.  This  means 
that  once  the  computer  determines  convergence  is  occurring, 
slow  convergence  is  increased  to  reduce  computer  time.  The 
fifth  and  final  entry  in  line  five  is  the  initial  temperature 
at  which  the  iterative  process  begins.  It  is  input  by  the 
user . 

Input  data  set  one  contains  temperature  dependent 
coefficients  and  is  not  used  in  this  model. 

Input  data  set  two  contains  up  to  50  constant  temperature 
inputs.  This  model  has  six  constant  ambient  temperatures  input 
by  the  user. 
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Input  data  set  three  involves  heaters  for  fast  warm  up  and 
is  not  used  in  this  model. 

Input  data  set  four  contains  all  pertinent  information 
concerning  the  ^i-node  equations.  Each  node  requires  two  lines 
of  data.  The  odd  numbered  lines  are  used  for  specifying  the 
nodes  that  interact  with  the  node  in  question  and  the  modes  by 
which  this  interaction  takes  place.  Consider  Figure  14,  line 
one  of  data  set  four  is  as  follows; 

6  7551  21  7521  111  7511  2411 

This  line  of  data  is  for  node  number  one  of  aspect  ratio 
selection  one.  The  first  entry  provides  the  number  of 
connections  to  that  node.  Entry  two  says  node  755,  an  ambient 
temperature,  is  connected  to  the  node  in  question  and  the  one 
indicates  that  the  connection  is  by  conduction.  The  rest  of 
the  entries  are  read  similarly  with  numbers  between  7511  and 
7561  indicating  a  conductive  connection  with  the  external 
environment.  External  heat  input  is  defined  by  entry  9991. 

The  even  lines  are  used  for  specifying  the  inter  node 
conductance  values  of  the  corresponding  entry  in  the  previous 
line. 

Input  data  set  five  is  used  only  if  there  are  unique 
exponents.  Therefore,  no  further  discussion  is  required. 

Input  data  set  six  contains  the  initial  temperature 
guesses  corresponding  to  the  number  of  nodes  receiving 
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secondary  heat  and  is  not  used  in  this  model  because  they  have 
been  set  by  the  last  entry  in  line  five. 

Input  data  set  seven  indicates  the  number  of  ..uu'.perature 
dependent  heat  input  curves  and  is  not  used  in  this  model. 

B.  FEATURES 

TMDL  is  presented  to  the  user  in  discrete  sections.  Upon 
entry  into  the  first  section,  the  user  is  offered  an  optional 
overview.  This  option  should  be  accepted  on  at  least  the  first 
run.  At  the  completion  of  the  overview  the  user  is  prompted 
for  a  data  file  name  and  title.  The  data  file  name  should  be 
changed  for  each  successive  run  because  TMDL  does  not  over 
write  existing  files. 

The  second  section  provides  prompts  for  the  physical 
characteristics  of  the  structure.  Initially,  the  user  is  given 
a  list  of  four  aspect  ratios  from  which  to  choose; 

1. )  10  by  24  nodes  provides  a  1:2.4  ratio  with  720  nodes. 

2. )  15  by  15  nodes  provides  a  1:1  ratio  with  675  nodes. 

3. )  12  by  20  nodes  provides  a  1:1.6  ratio  with  720  nodes. 

4. )  8  by  30  nodes  provides  a  1:3.75  ratio  with  720  nodes. 

It  should  be  noted  that  up  to  30  additional  nodes  may  be  added 
depending  on  the  existence  of  the  mounting  ear.  Figure  16 
displays  the  specific  geometry  with  definitions. 

After  selecting  the  desired  aspect  ratio,  the  user  must 
specify  the  use  of  either  SI  or  English  units.  The  input  of 
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BOTTOM 


Figure  16 


data  must  remain  consistent  with  this  selection.  At  this  point 
the  input  of  structure  characteristics  are  requested.  These 
entries  consist  of  the  length,  width,  thickness,  and  thermal 
conductivity  of  each  layer.  There  is  ample  opportunity  for 
correction  or  alteration  throughout  TMDL. 

Section  three  requests  the  input  of  initial  and  ambient 
temperatures.  TMDL  has  six  ambient  temperatures. 

The  final  section  allows  for  the  input  of  external  heat 
sources.  TMDL  provides  for  heat  input  into  the  upper  substrate 
surface  only.  There  are  four  methods  of  external  heat  input 
from  which  to  choose.  The  first  alternative  allows  for  a  total 
rate  of  heat  applied  to  the  upper  surface.  An  entry  for  this 
alternative  would  be  divided  by  the  nodes  and  distributed 
appropriately.  Option  two  provides  for  the  entry  of  average 
heat  per  unit  area.  The  third  alternative  givej.  the  user  the 
freedom  to  enter  heat  in  specifically  selected  nodes. 

Option  four  offers  no  heat  input. 

At  the  conclusion  of  these  entries  an  output  data  file 
will  be  created  in  the  proper  format  and  placed  in  a  file 
named  by  the  user. 
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APPENDIX  A 


TMDL  LISTING 


S LARGE 
C 

C  TITLE;  MODEL  BUILDER 
C  AUTHOR:  LT  PATRIC  ROESCH 
C  DATE:  15  MARCH  1991 

C  COMPILER;  MICROSOFL  VERSION  4.01 
C  LINKER:  MICROSOFT  VERSION  3.5.5 
C 

C  DEFINE  REAL  VARIABLES 

REAl.  SL.SW,DELS.KS,EW,EL,DELL.KL,CL,CW,DELC.KC,IT,UPR  r.LWRT,LT,RT 
+  .FT,BT,THEAT,THPN,AHEAl,NHEAl.DELX,DELY.SYLR.SXrB,SXX,SYY,SZE,SZT 
REAL  EXEB,EYY.EXX.EZC,L/\R,CYLRE.CYYE.CXFBE,CXXE,CZEAR,CYLR,CYY,CXX 
+  ,CZB,ACC.DAMP,CONF.\C.EYLR 
C 
t: 

C  DEFINE  IN  TEGER  VARIABLES 

INTEGER  NPL,NWIDE.NDEEP,NUM.NUMA,CH,H,TOTNOD,NN.NC,N,IB,IE,IC,ID, 

+  COUNl,CONTEMP.USEL,ZER,NMAX.TMAX,HTRS,DLD2.D3.D,5.D4,D6,D7,MAXIT 
C 

C  DEFINE  /\I.L  ONE  CHARACTER  VARIABLES 

CHARACTER*  1  ANS,ANSN,ANSA,ANSSL.ANSSW,ANSST.ANSSK,ANSE, ANSEL. 
+.ANSEW,ANSET.ANSEK.ANSC,ANSCL.ANSCW.ANSCr./\NSCK,ANSr.ANSTLANSTU, 
■h.VNS  LLR.ANS'I  L.ANS  TR.AN.'^  I  K.-YNm  rB.ANSH,ATH,ANSHA.,\HNALLECr.SELH 
+  ,ANSWER.DLLr 
C 

C  CHARAC'lER  VARIABLES  OF  MORE,  l  llAN  ONE  POSI  IION 

CHARACTER  DATAF*70,UK*1LUT*LUH*13,UAH*6,NAME*6,UL*2 
C 

C  DEFINE  MATRICES 

REAL  HEAT(30,1.5),COEF(740,9) 

C 

INTEGER  IH(240),JH(240),NCON(740,9) 

C 

c: 

C 

C  PROVIDE  THE  USER  WITH  A  PROGRAM  OVERVIEW 
( 

HKK)  CALL  CLS 
WRITE!’, 1001) 

1001  FORMAT!///.' 

+  I  HLS  PROGRAM  WAS  WRLITEN  Id  INTEGRA!  E  W1  HI  I  XIS  I  ING  ,/. 

+  I  HERMAL  ANALYSIS  SOFTWARE  AND  TO  REDUCE  THE  AMOUNT  ,/, 

+  OF  [IMF  REOUIRED  FOR  DATA  ENTRY. 

+  WOULD  YOU  LIKE  AN  OVERVIEW  OF  THE  PROGRAM  PRIOR  TO  ',/. 

+  BEGINNING.’  ENTER  Y  FOR  YES  OR  N  FOR  NO.  ■.2X,\) 

READ!*. ’002)  OVR 
too:  FORMA  I  (A1) 
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IF  (OVR.EQ.’Y'.OR.OVR.EQ.’y’)  THEN 
CM.L  CLS 
\VRirE(M003) 

1003  F0RMAT(  «..***»»««OVERVlEW‘*‘******  ■,//, 

+  rms  PROGRAM  PERFORMS  A  NODAL  /VNALYSIS  OF  A  THREE  V, 

+  LAM  RED  CONDUCnVE  S  lRUCrURF  (SUBS  TRATE.FPOXY.AND  . 

+  C' ARRIFR).  I'HE  OU  I'PU'l'  CONSISTS  OF  UP  TO  740  COEl  l  l-  ,  . 

+  c  if;n  IS  iHAT  CON iRiBU  n.  ro  nil  oi n  rminaiion  oi 

+  IHF  rEMPERAI  URE  DISTRIBLTION  OF  A  MICROCIIU  U1 1 
+  WHEN  FED  IN  TO  PHE  LHERM/U.  ANAl.Y/F.R.  „  . 

+  rilE  FOLLOWING  IS  /\N  ORDl-RED  OUTLINE  OF  THE  MAJOR 
+  SECTIONS  OF  THIS  PROGRAM  AND  WHAT  ENTRIES  ARE 
+  REQUIRFi'  OF  THE  USER. 

+  NO  I  L  I  HA  L  ALL  EN  TRIES  WILL  BE  IN  UPPER  CASE  LE  TTERS. 

+  A.)  DATA  FILE 

+  1.)  ENTER  THE  TITLE  OT  THE  OUTPU  T  DATA  FILE.  ,/, 

+  2.)  ENTER  THE  DATA  EILE  NAME.  THIS  PROGRAM  WILL  NOT’,/, 

+  ERASE  OR  WRITE  OVER  EXISTING  EILES.  'Jl, 

+  PLEASE  PRESS  ENTER  WHEN  READY  TO  CONTINUE.  ■.2X,\) 

READ(  M()02)  ANSWER 

c: 

CALL  CLS 
WRTTE(M004) 

1004  FORMAT(/V,' 

+  B.j  STRUCTURE  PHYSICAL  CHARACTERISnCS 
+  1.)  SELECT  UNIT  TYPE  (SI  OR  ENGLISH) 

+  2.)  SELECT  FROM  FOUR  ALTERNATIVES  THE  DESIRED  NODAL’,/, 

+  ASPECT  RATIO.  ’,/, 

+  •  3.)  ENTER  EACH  LAYER  CHARACTERIS'nCS.  SUBSTRATE  IS  ’,/, 

+  ■  lilE  TOP  LAYER,  FOLLOWED  BY  EPOXY  AND  THE  CARRIER',,'. 

+  LAYERS.  RESPECTIVELY.  nJE  CARRIER  LAYER  MAY 

+  CON  TEMN  .'\N  EAR  ON  THE  FRONT  /\ND  BACK  SURFACES 

+  IF  SPECIFICATIONS  REQUIRE  IT.’,//. 

+  PI.EASE  PRESS  ENTER  WHEN  READY  TO  CONTINUE. 
RF.AD(M002)ANSWF.R 
C 

CALL  CLS 
WRTTE(MOO.S) 

100.>  I-ORMAT(,///. 

+  C.)  INITIAL  AND  AMBIENT  TEMPERATURES  ’,/, 

+  1.)  ENTER  INITIAL  CHIP  AND  CHASIS  TEMPERATURE.  ’,/, 

+  ASSUME  STEADY  STATE  FOR  INITIAL  TEMPERATURE.’,/, 

+  2.)  ENTER  AMBIENT  TEMPERATURES  FROM  ALL  SIDES  OF  THE’,/, 

+  '  STRUCTURE.’.//, 

+  ■  D.)  HEAT  INPUT  ’,/, 

+  ■  HEAT  INJECTION  OCCURS  ONLY  ON  THE  UPPER  SUREACE  OF  ’,/, 

+  THE  SUBSTRATE  LAYER.  THIS  PROGRAM  SUPPLIES  THE  USER  ',/, 

+  FOR  ALTERNATIVE  METHODS  FOR  ENTERING  HEAT;  ,/, 

+  1.)  TOTAL  HEAT  OVER  SUFACE  ’,/, 

+  2.)  AVERAGE  HEAT  PER  UNIT  AREA  ’,/, 

+  3.)  INPUT  HEAT  NODE  BY  NODE  ’./, 

+  4.)  NO  HEAT  INPUT  ’,//, 

+  THIS  COMPLETES  THE  PROGRAM  OVERVIEW.  PLEASE  PRESS  ENTER  ’,/, 
+  TO  CONTINUE.  ■,2X,\) 
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READ(M002)ANSWER 

C 

ELSEIF  (OVR.EQ.’N’.OR.OVR.EQ.  n’)  THEN 
GOTO  1007 
ELSE 

GOTO  1006 
1007  ENDIF 

c: 

C 

c: 

G  INITIAI.IZE  MATRICES 
C 

DA'I'A  HEAL  M.SO’O.O/ 

DATA  IH  /240*0/ 

DATA  JH  /240*0/ 

(' 

G  VVMilAULE,  GONSI  ANT.  AND  SI  RING  DEEINHION 
(  PHYSIGAL  GHAIMCrr  RISriGS 

G  SL,EL,CL  -  SUBSTRATE.EFOXY,  /\ND  GARRIER  LENGTHS 

(  S\V,EVV,C\V  -  SUHS  I  RATE.EPOXY.  AND  GARRII  R  WIDTHS 

G  KS.KE.KC  -  SUBSRA'PE.EPOXY.AND  GARRIER  Till  RMAL  CONDUGTIVITIES 

G  DELS.DELE.DELG  -  SUBS  I  RA  LE.EPOXY  AND  GARRIER  THIGKNESS 

C  DELX  -  SL/NDEEP 

(  DELY  -  SW/NWIDE 

C  NPL  -  NUMBER  OF  NODES  PER  A  LAYER 

G  NWIDE  -  NUMBER  OF  NODES  WIDE 

C  NDEEP  -  NUMBER  OF  NODES  DEEP 

C  UL  -  UNITS  OF  LENGTH  (SI  OR  ENGLISH) 

C  UK  -  UNITS  OF  THERMAL  CONDUCTIVITY  (SI  OR  ENGLISH) 

C  EAR  -  DEPTH  OF  THE  EAR  (SL-CL)/2 

C 

G  INITIAL  AND  AMBIENT  TEMPERATURES 

C  IT  -  INITIAL  CHASIS  AND  CHIP  TEMPERATURE 

G  LT  -  LEFT  SIDE  AMBIENT  TEMPERATURE 

G  RT  -  RIGHT  SIDE  AMBIENT  TEMPERATURE 

t:  FT  -  FRONT  AMBIEN  T  TEMPERATURE 

G  BT  -  BACK  AMBIEN  T  TEMPERA  TURE 

G  UPRT  -  UPPER  AMBIENT  TEMPERATURE 

G  LWRT  -  LOWER  AMBIENT  TEMPERATURE 

(  U  T  -  UNI  TS  OF  TEMPERATURE  (GENTIGRADE  OR  FARENHIE  T) 

G 

G  HEAT  INPUT 

(  THEAT  -  TO  TAL  INJECTED  HEAT 

(  THPN  -  TOTAL  HEAT  PER  NODE 

(  AHEAT  -  AVERAGE  HEAT  OVER  A  GIVEN  SURFACE 

(  NHEAT  -  HEAT  PER  NODE  INJECTED  NODE  BY  NODE 

(  UH  AND  UAH  -  UNITS  OF  HEAT  (SI  OR  ENGLISH) 

(  NUM,NUMA.CH,H  -  DUMMY  VARIABLES  USED  TO  CREATE  THE  VECTORS 

(  IH,  &  JH  WHICH  ARE  USED  lO  RELA  TE  NODE 

(  NUMBER  TO  MATRIX  POSTTION 

(  TOTNOD.NN.NCN  -  VARIABLES  USED  TO  ALLOW  HEAT  INPUT  NODALLY 

C:  IFI.JH  -  VECTORS  USED  TO  CORRELATE  NODE  NUMBER  WITH 

G  MATRIX  POSITION 

G  HF;AT  -  MATRIX  USED  TO  CONTAIN  HEAT  INPUTS 

( 
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COEFFICIENT  DlTlNtTIONS  Y  IMFiJF:S  WIDTH  &  X  IMPI.IFS  DEPTH 
Z  IMPLIES  HIEGH  T 

S  IMPLIES  SUBSTRATE,  E  IMPLIFTS  EPOXY 
C  IMPLIES  CARRIER 

SYLR  -  PROVIDES  COEFEICII  N  T  LOR  LEFl'  OR  RICH  T  EDGE  NODES 
TO  THE  EXTERNAL  NODE.  SYLR  IMPLIES  SUBSTRATE, 

EYLR -  EPOXY 

CYLR  -  CARRIER/  NO  EAR 

CYLRE  -  CARRII:R  W/  EAR 

SXFB  -  PROVIDES  COEFFICIENT  FOR  FRONT  AND  BACK  EDGE.  NODES 
TO  THE  EXTERNAL  NODE. 

EXFB -  EPOXY 

CXFB  -  CARRIER/  NO  EAR 

CXFBE  -  CARRIER  W/  EAR 

SYY  -  PROVIDES  INTERNAL  COEFFICIENT  IN  THE  ’Y'  DIRECTION 
EYY -  EPOXY 

CYY  -  CARRIER  W/NO  EAR 
CYYE  -  CARRIER  W/  EAR 

SXX  -  provide;  internal  COI.I  I  ICIENT  IN  THE  X'  DIRECTION 
EXX -  EPOXY 

CXX  -  CARRIER  WNO  EAR 
C.XXE  -  C/\RRIFR  W'  FAR 


SZT  -  COEFFICIENT  FOR  SUBSTR.ME  TO  UPPER  EX  TERNAL  NODE 

SZE  -  COEFFICIENT  FOR  SUBSTRATE  TO  EPOXY 

EZC  -  COEFFICIENT  FOR  EPOXY  TO  CARRIER 

CZEAR  -  COEFFICIENT  FOR  EAR  TO  LOWER  EXTERNAL  NODE 

CZB  -  COEFFICIENT  FOR  CARRIER  TO  LOWER  EXTERNAL  NODE 

CHARACTER 

ALL  ONE  CHARACTER  STRINGS  EXCEPT  ’UT  ARE  SIMPLE 
YES.  NO,  OR  SELECTION  ANSWERS  AND  DO  NOT  REQUIRE 
EXPLANATION. 

DATA  FILE 

DATAF  -  TITLE  FOR  DATA  FILE,  LINE  ONE 
NAME  -  DATA  FILE  NAME 

IB,IE,ID,IC  -  COUNTERS  FOR  HLLING  DATA  FILE 
COUNT  -  TOTAL  NUMBER  OF  NODES 

ACC  -  MINIMUM  ERROR  CRITERIA  AT  WHICH  CALCULATIONS  CEASE 
DAMP  -  DAMPING  FACTOR  FOR  PREVENnON  OF  TEMPERATURE 
OSCILLATIONS 

CONFAC  -  CONVERGENCE  FACTOR  USED  TO  INCRIiASI  DAMPING 
CONTEMP  -  TMDL  HAS  6  CONSTANT  /XMBII  N  T  TEMPI  RA TURFS 
USEL  -  TELLS  ANAI.YZER  WHAT  1 YPF  OF  UNI  TS  ARI-  IN  USE 
ZER  -  DUMMY  VARIABLE 

NMAX  -  MAXIMUM  NODES  PROGRAM  IS  (  APABLl  OF  EVALUATING 
I  MAX  -  M/\XIMUM  NUMBER  OF  CONS  IAN  T  TEMPERA  TURES  PROGR,U1 
C/\N  ACCEPT 

H  TRS  -  NO  T  APPLICABLE  TO  THIS  PROGRAM  -  NO  EXPLANATION 
NECCESSARY 
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('  D1-"  -  data  sets  in  use; 

(  NUVXIT  -  M^VXiMUM  ITERAllONS  COMPU  TlTi  WILE  IT-RIORM  TO 

(  ACHIEVE,  SOEU  nON 

c 

(  HE(ilN  I’ROORAM 
( 


(  .\EE  C  I  S 


\VRllEi',40«) 
4()(.  1C3R.\E\1(  , 


-h  ■ 
-t-  ■ 
+  ' 


Pi, EASE  Sl  l.l-CI  UPPER  (  ASE  1,1  El!  RS 
PRIOR  K)  HEOINNINC). 


( 


*««*«*««««*•««•*«**«••««««««*«*«««  «*W>ti****«*)k*«»1‘*‘ 


WRirE(*.2()l) 

2i)l  t\)rmai{  ,'  ■,/, 

+  PRIOR  TO  EN'n;RING  DA  I  A  INTO  THIS  PROGRAM  ENSURE  THAT’,/, 

+  Y(3U  HAVE  A  DRAWING  OF  YOUR  DESIGN  AND  ALL  PERTINENT  V, 

+  DATA.  './■/, 

+  PRESS  EN  rE;R  WHEN  RliADY  lO  CONTINUE.',2X,\) 

re;.\D(*,304)answe:r 

(■ 

(  ALLOW  USE.R  K)  N/VME  THE  DAI  A  I  ILE 
C  ALL  CLS 

(■ 

WRriL,(’.'()l) 

'1)1  IOR.MAK  ,  nils  PROGRAM  CREAH  S  AN  OU  ITU  I  DAl  A  FILE  I'OR  EN  TRY 
+  INTO  nil  ',. 

4-  E.XLSriNG  I'HERM/\J.  ANALY/I  R.  FURIHERMORE.  THIS  PRCXiRAM  DOHS  ,. 
4-  NOT  ERASE  OR  WRITE  OVER  THE  EXISTING  DATA  FILE,  THEREFORE'.. 

4-  nil  USER  WILL  name;  Tin  D.VTA  file  for  eac  h  run  of  THIS  ... 

4-  PROGRAM,  the:  TILE  NAME  IS  LIMI  TED  TO  SIX  Cl lARAC'TFRS, 

4-  PLEASE  ENTER  THE  DI  MRI  D  DATA  FITE  N.AMi-.  ■.2X.  ) 

READ  '02.  N/^^IE 
't)2  EORMATlAO) 

( 

( 

( 

(  ALLOW  USER  TO  PROVIDE  TTri  E  LINE  TOR  DATA  FILE 
CALL  CLS 
WRITE!  *,202) 

202  FORMATE/'/.'  ENTER  THE  DESIRED  mT.I-,  TO  HE  PI.AC/ED  ON  LINE 
4-  NUMBER  ONE  OF  THE  OUTPUT  DATA  ETLE.'... 

+  ■  ',2X,\) 

RE;/\D  l.DATAF 
1  FORMAT(A70) 

( 


C  ALL  CLS 


( 

( 

{  M  PPLY  A  LIS  TINC)  (JE  AC  C'l  1’ TABLE;  ASPEC  T  RATIOS 
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C' 

\vKrii;c,203) 

:u3  10RMAT(  .  rUH  mLl.(3WING  ASPECT  RAIIOS  ARP,  AVAILABLE.  ) 

WRI  rEC,204) 

2lt4  {ORMATi  1.)  10  BY  24  PROVIDF-S  A  1:2.4  RATIO  Wi  l  l  I  740  NODES. 

+  .  2.)  I.^  BY  IS  PROVIDES  A  1:1  R.ATIO  WITH  675  NODES.’,'. 

+  .7.)  12  BY  20  PROVIDES  A  1:1.6  RATIO  WITH  720  NODES.’., 

+  .  4.)  8  BY  30  PR0V1DL;S  A  1:3.75  RATIO  WITH  720  NODES.’, 

+  .  PLEASE  SELECT  A  NUMBER  ONE  LHROUGH  FOUR.  .2X,\) 
READ(L.3020)SELE(T 
3020  F’ORMAKAij 
C 

(’  USER  CAN  CON  llNUE  OR  MARE  ANOlllLR  SEI  EC’llON 
(  ALL  CLS 

IP  (SI  l.EC’l.NE  1  AND.Sl  LI  (  I.NE.  2  AND.SELE.Cl.Nl  .’3’.AND.SELECr. 

+  Ni;.’4’)  niEN 
CjO  PO  7 
ELSE 

111  WRI  I  E{*,3)  SEI.l  ("I 

EORMAK  YOU  SEI  I  (”IED  NU.MBI:R ’.Al. OF  niE  EOLLO 

+  VMNCi  4  AL  ITRNAI  IVES.  „  . 

+  1.)  10  BY  24  EOR  A  1:2.4  R.-VEIO  ,  . 

+  ■  2.)  15  BY  15  EOR  A  1:1  RAIIO’.,, 

+  .7.)  12  BY  20  E'OR  A  1:1.6  R.\TIO’.,. 

1-  4.)  8  BY  .30  EOR  A  1:.3.75  R.-Vl  10„  ) 

C 

ENDIE 

C 

5  WRITE!  *.205) 

2(\5  I ORMAL!’  IS  THIS  THE  DESIRED  SEI.ECTION.  ENTER  Y  EOR  YES  AN 

+  D  N  EOR  NO.  ’,2,\,  ) 

READ(*.  8)ANS 
8  EOR.MAT(Al) 

(’ 

C 

IE  (ANS.EQ.’N  )  THEN 
GO  EO  7 

ELSLIE  (ANS.EQ.’Y’i  Elll  N 
GO  l’O  0 
ELSE 

(  ALL  (  LS 
(jOEO  10 
M  ENDIE 
C 
c 

CALL  CLS 

( 

C 

(  ,\E*I  ER  CONEIRM.A  nON  OE  SELECTION  FILL  C(3NST/VNTS  WITH  APPROPRIATE 

(  V..\LUES 

IE  (SELEC’I.EO.  I  )  THEN 
NWIDE  =  10 
NDEEP  =  24 
NPL  =  240 
C 
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I'l.SlIlF  (SELlX  l  .l-.U.  2  )  I  IILN 
NWIDE  =  15 
NDtbP  =  15 
NPL,  =  225 

( 

1  LSEII-  (SELhCl.l  U.'5  )  lllEN 
NWIDE  =  12 
NDEEP  =  20 
NPL  =  240 
C 

El.SElF  (SELECEF.O.'4  )  THEN 
NWIDE  =  8 
NDEEP  =  30 
NPL  =  240 
ENDIF 

c: 

C  MAKE  DESIRED  UNFL  Sf-LECTION 
C 

2  WRrrE(*,206) 

206  FORMAK  "  .  PHLS  PROGRAM  IS  CAPABLE  OF  OPERATIONS  IN  El 

4- 1  HER  SI  OR 

+  ENGLISH  LNl  I  S.  AFTER  THE  SI  LI  (TION  OF  FHE  UNI  IS,  ALL  ■,/. 

4-  EINTRIL.S  MUS  I  BE  C  OMPAHBLi:.  PLEASE  MAKE  YOUR  SELECI  ION.',/ 

4-  S  FOR  SI  NO  FATION', , 

4--  L  FOR  FN(. FISH  NOI  AIION  ■.2X.,) 

RE.AD|*,302)ANSN 
,502  F(3RMAT(A1) 

(■ 

(■  CHECK  FOR  CORRLCl  UNFl  SLLLCHON 
4  IF' (/\NSN. EXE  S' )  FHEIN 

WRITE!  *,207) 

20’  FORMAT!///,.'  YOU  SELECFED  SI  NOTATION.') 

ELSEIF  !ANSN.EQ.'I-')  THEN 
WRITE!  *.208) 

208  FORMAT!/// .'  YOU  HAVE  SELECTED  ENGLISH  NOTATION.  ) 

ELSE 
CALL  CLS 

GO  ro  2 

ENDIF 

(' 

W'Rri'E!’.209) 

20^  FORMA  T!  18  THIS  THE  DESIRED  SELECTION  .’  EN  TER  Y  FOR  YES  AND 
4-N  FOR  NO.  '.2X./) 

RFAD!*,303)ANSA 
,503  FORMAT!  Al) 

( 

( 

(  SHOULD  I  CONTINUE  OR  RESF;I  FCT  UNI  TS 
IF'  (ANSA.EO  'Y  )  THEN 
GO  TO  6 

)  I  SFIF  I  ANSA.FO  N  1  TH)  N 
(  ALL  CLS 
GO  TO  2 
ELSE 

CALL  CLS 
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WRITE(*,210) 

:U)  FORMAT(;/.,  ,  ,  THIS  PROORAM  IS  CAPAlil.h  OIOPF.RA  I  IONS  IN 
+  FITHHR  SI  OR 

+  ■  FNGl.ISll  L'M  IS.  Al'l  l.R  lllF  SF:i.F-C:riON  OF'  UNFl'S.  ALL',/. 

+  LN  I  RII  S  MUS  I  BL  COMP.A  I'lBLI-.  IT,I;aSI-;  MARI-  YOUR  SFLFC'  I  ION.', 
-r  .  S  FOR  SI  NO  r, VI  ION.' 

+  F  FOR  I  NOLISH  NOI  AIKJN.') 

(' 

GO  I O  4 
(.  I  N  OIF 
(■ 

(  ALL  CT.S 


U 

^  :M  «  V  >!•  W  It  «  «  «  «  «  «  «  «  «  S  --V  --S  1:  s  <5  i:  -1;  >  -V  Is  i: 

LN  I  I  R  (  HIP  CHARACTHRISTICS  FOR  EACH  LAYER**** . * 

(' 

C 

l^i^yER* 

C 

wRrrE(*,2ii) 

211  FORMAT! 


SUBSTRATE  CHARACTERISTICS** 


C 

C  PROVIDE  CORRFC  l  UNFl'  ABBREVIATIONS 
IF  (.ANSN.F.IT  S'i  IHF.N 
WRITE! '.212  i 

212  FORMA IV  ALL  I  N  FRIES  ARE  IN  SI  NOTATION.',,) 

UL  =  'em' 

UK  =  'Waits  tm.C  ' 

U'F  =  '(" 

I  I.SFIF  (,\NSN.l  O.  I  ')  Fill  N 
WRFFE(*.21.^) 

21.^  FORMAK  ALL  EN  FRIES  ARE  IN  ENGLISH  NO  rAnON.  .,0 
UL  =  'in' 

UK  =  BiuTir  lal  ' 

U'F  =  I" 

ENDIF 

(' 

C 

(' 

WRITE(*,220)UL 

220  FORMATE  ENTER  SUBSTRATE  LENGTH  (’.A2.');  .2X,\) 

RE^VD  ‘.SL 

(' 

WRFFF(*,216)UL 

210  FORMAl'f  ,'  ENTER  SUBSTRATE  WIDTH  r.A2.');  ■,2X,\) 

Rf  AD  ‘  SW 

( 

WRI  Fit*  2r  )UL 

:r  FORMAL  I  NiFR  .oUBSFRAFF,  FHICKNESS  t'./\2.'):  ■.2X.  ) 
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READ  ’.DI-l.S 


WRnE(*.21S)UK. 

21.S  f-ORMAT(/,'  ENTER  SURSfRAlE  rilE.RMAL  CONDUd  lVlTY  C.An.'):  ’,2X,\) 
READ  EKS 
C 

31  CAEE  CES 
C 

C;  M/VKE  ANY  C  l  lANGES  C3R  CIORRECTIONS  TO  SURSRATE  EN  TRIES 


214 

20 


C' 


WR[TE(E214) 

EORMAl  l  ,  YOU  HAVE  MADE  CHE  I'OELOWING  ENTRIES  FOR  THE  SUB 
+  STRATE. 

ERIN  I  N'  1.)  l.ENG  I'H  '.SE 
PRINT*,'  2.)  WIDTH  '.SW 

PRINT  E  3.)  IKK  RM  -SS  .DEES 

PRINT*.'  4.)k  ',KS 


WRTTE{*.2ES) 

21.3  FORMAT(;.  DO  YOU  WISH  TO  MARE  ANY  CHANGES?  SEEECrT  Y  FOR  YES 
4-  AND  N  FOR  NO.  ■.2\,o 
REAI3(*..W)ANSS 
.304  EORMAKAl) 


IE  (ANSS.EO.'Y')  THFIN 
12  CAEE  CES 
WRITE!*, 219) 

219  FORMAT!/';, 0 

PRIN  T  *,  THE  CURREN  T  ENTRY  FOR  LENGTH  IS  '.SL.'  '.UL 
WRITE!  *.2190) 

2190  FORMAT!/,'  WOUED  YOU  EIRE  TO  CHANGE  THE  EENG  TH?  (Y  OR  N) 

+  ■.2X,o 

READ(*,.304)ANSSE 
PRINT  * 

IF  !ANSSL.F0.'Y')  then 
WRTTE!*,22I)UE 

221  FORMAT!  '  EN  TER  SUBSTRATE  EENG  TH  (■./\2.'):  ■.2X.,) 

READ  ‘.SI 

EESFIF  I.ANSSE.EQ.  N  )  Till  N 
GOTO  14 
FESF 
CjOTO  12 
ENDIF 

( 

14  CAEE  CES 

WRITEC.219E) 

2191  FORMAT!,' ,'//■) 

PRINT  *.'  THE  CURRENT  ENTRY  FOR  WIDTH  IS  '.SW.'  '.UE 
WRITE!  *,222) 

222  FORMAT!  ,'  WOUED  YOU  EIRE  TO  CHANGE  THE  WIDTH?  !Y  OR  N)  ' 
+  2X..) 

READ! *.2 192)  ANSSW 

2192  F(3RMAT!A1) 

IF  !ANSSW.E0.  Y  )  THEN 
WRITE!*.223)UE 
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IM  I  R  SL  liS  1  RAl  i;  WID  I  H  (  .a:.  ):  ■.2X,,) 


::2  tORMA'K  . 

RHAD  ‘.SW 
Rl.SMR  (ANSSW.I  CX'N  )  THI  N 
(ioro  12 

l-LSI- 
GOTO  14 
f-NDlI- 
G 

c: 

13  CALL  CI.S 

WRITE!  *,2193) 

2193  FORMAT) ///O 

PRINT  THE  CURRENT  ENTRY  FOR  THICKNESS  IS  ’.DELS,’  ’.UL 
WRITE) ’,224) 

224  FORMAT)/,’  WOULD  YOU  LIKE  TO  CHANGE  THE  THICKNESS?  (Y  OR  N) 

+  '.2X,\) 

READ)*, 2194)  ANSST 

2194  FORMATIAI) 

IF  )ANSST.EO.  Y  )  THEN 
WRni)*,22.3)UI 

223  lORMAT).  I  N  PER  SUIISTRATE  I  HICKNESS  (G\2,  ):  ■.2X. 

+  ) 

Ri:,\D  *.DELS 

ELSEIt-  (ANSSTLO.  N  i  THEN 
GOrO  13 
ELSE 
GOTO  13 
ENDIF 
C 
C 

13  CALL  CLS 

WRITE)*,219.3) 

2193  FORMAT)///) 

PRINT  THE  CURRENT  ENTRY  FOR  THERMAL  CONDUCTIVITY  IS  ’,KS 
+  .■  ’.UK 

WRITE)  *.226) 

226  FORMAT)/,’  WOULD  YOU  LIKE  TO  CHANGE  THE  THERMAL  CONDUCTIVIT 
+  Y?)YORN)  ’,2X.\) 

READ)*, 2196)  ANSSK 
2196  FORMAT)Al) 

IF  )ANSSK.EQ.’Y  )  THEN 
WR1TE)’,227)UK 

22'  FORMAT)  ,'  EN  PER  SUHS  PRAPE  THERM/\1  CONDUC’ PIVPPY  )',A1 

>2.'):  ■.2X.\) 

READ  ’.KS 

ELSEIF  lANSSK.EO.  N  )  IHP.N 
GOTO  16 
ELSE 
GOTO  13 
ENDIF 
C 

(  ALLOW  ANO'PHER  REVIEW  OF  SUBSPRAPE  DATA  EN'PRIES 

C 

16  CALL  CLS 
WRITE)  *.2281 
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22cS  FOI^MATf///;  YOU  MAVI-.  MADE  TIIE  I'OLEOWING  CORRECTIONS  TO  1' 

+  HE  SUBSTRATE  EN  TRIES.  /) 

CO  TO  20 
C 
( 

t.l.SEll-  (ANSS.ECJ.  N  )  Till  N 
(lOK)  r 
TEST 
GO  TO  .M 
ENDIF 
C 

c 

characteristics************ . . 

c 

c 

17  call  CLS 
WRITE!  *,229) 

229  FORMATE//,' . * . . . . . 

+  ' . . CHARACTERISTICS . 

^♦♦■f****’  j 

+  ■  NOTE  THAT  SUBSTRA  TE  AND  EPOXY  LENGTH  AND  WIDTH  ARE  EQUAL.  ',/ 


+  ) 

C 

WRTTE(*,230)UL 

230  FORMAT!/.'  EN  TER  EPOXY  THICKNESS  !'.A2.'):  '.2X.\) 

READ  *.DELE 

WRTTE!*.23I)UK 

231  FORMAT!  .  L.N  1 1  R  I  POX'!  IIII  R.MAI  CONDUCTIVITY  !'..\12.  ):  ■.2X.,) 
READ  ’.KE 


C 

EL  =  SI. 

IVV  =  SW 
C 

C'  REVIEW  EPOXY  EN  TRIES 
(' 


30  CALL  CLS 


WRITE! ’,232) 

232  FORMA  T!  7,  '  YOU  HAVE  MADE  THE  FOLLOWING  ENTRIES  FOR  THE  EPOXY 


+  LAYER.'./') 

32  PRINT  *.'  1.)  LENGTH 

PRINT  2.)  WIDTH 
PRINT*.'  3.)'raiCKNESS 
PRINT  *,’  4.)  k 


,EL 

'.EW 

DELE 

'.KF 


(' 


WRITE!  *.233) 

233  FORMAT!  .'  NOTE  THAT  CHANGING  EPOXY  LENG  TH  OR  WIDTH  AI.SO  CHANG 
+  L;S'.,.'  SUBSTRATE  LENGTH  OR  WIDTH.',,/. 

+  DO  YOU  WISH  LO  MAKE  ANY  CHANGES  TO  TFIE  EPOXY  ENTRIES.’  !Y  OR 


+  N)  '.2X.,) 

READ!  *.219^1  ANSL 
219^  FORMATlAl) 


C 
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C'  MAKE  ANY  CHANGES  OR  CORRECTIONS  TO  EPOXY  ENTRIES 
C 

IE  (ANSE.EO.  Y  )  THEN 
::  CALI.  CLS 

WRITE!*, 2198) 

2198  EORMAT(;'.0 

PRIN  T  *.’  THE  CURREN  T  EN  TRY  EOR  LENC.1II  IS  '.EL.'  '.UL 
WRITE)  *.234) 

234  TORMAT)  ,  WOULD  \  OU  LIKE  TO  CHANGE  THE  LENGTH?  (Y  OR  N) 

+  ',2X,  ) 

READ)*, 2199)  ANSEL 
21W  EORMAT(AI) 

C 

IT'  (/XNSEI  .EO.'Y  )  THE.N 
WRTTL(*,23.3)UL 

233  E(3RMAT(,,'  1  N  il  R  epoxy  LI  NGTH  (•.A2.-);  •,2X,,) 

READ  Mil. 

SL  =  EL 

EILSEIF  (ANSLL.LO.  N')  THEN 
GO  TO  2 1 
ELSE 
GCTTO  22 
21  ENDIE 

C 

24  call  CLS 
WRITE!*, 6G1) 

901  EORMAT!////) 

PRINT  *,•  THE  CURRENT  ENTRY  FOR  WIDTH  IS  '.EW,'  ’,UL 
WRITE!  *,236) 

236  F-'ORMAT!-,'  WOULD  YOU  LIKE  TO  CHANGE  THE  WIDTH?  !Y  OR  N) 

+  '.2X.'\) 

READ!*.602)ANSEW 

602  FORMAT!  Al) 

C 

IE  (/VNSEAV.EQ.  Y  )  THEN 
WRITE)  *.23 1’jUL 

23*  EORMAT)  ,  EN  TER  EPOXY  WIDTH  )  ,A2.  );  .2X.  ) 

read  *,EW 
sw  =  e;w 

ELSE'IE  lANSEW.EU.  N  )  THEN 
GO' TO  23 
ELSE 
GOTO  24 
23  ENDIF 

C 

26  C'ALL  CLS 
WRITE!  *,603) 

603  FORMAT!/;,) 

PRINT  *,•  THE  CURRENT  ENTRY  FOR  THICKNESS  IS  .DELE,’  ’.UL 
WRITE!  *  238) 

238  EORMAT!  .  WOULD  YOU  LIICE  TO  CHANGE  THE  THICKNESS?  !Y  OR  N) 
4  ■,2X,,) 

READ!*,604).ANSET 

604  F0RMAT!A1) 

(' 
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IF  (,\NSFT.FU.  Y  )  IIIFN 
WRn'F:(M>F)L'l 

l  ORMAi  f  .'  1  N  Ti  R  I  POXY  niK'KNi  ss  (  .a:.  );  ■.:x.  ) 

RFAH  M)l  l  F 

Fi.sFii-  (ANSF  r.i':u.'N  )  tiii-;n 

C.O  TO  25 
ll.SF 

cio  rc)  20 

25  FNDIF 

(' 

2S  CALL  CI  S 
\vRrrE(*,w)5) 

(>05  FORMA  r(w,) 

PRINT  Y'  THE  CURRENT  EN  TRY  FOR  THERMAL  CONDUCTIVI  TY  IS  MCE 
+;  '.UK 

\VRTTE(*.240) 

241)  FORMA'T(;,'  WOULD  YOU  LIKE  TO  CIEYNGE  THE  THERMAL  CONDUCTIVTT 
+  Y.’(YORN)  'MX.) 

READ!  YW)6) ANSF  K 
l-ORMATiAl) 

( 

IF  (ANSFK.I  0.0  )  THI  N 
WRTTE(Y241)UK 

241  FORMAK,  FNFFRFPOXY  1 1 IFRMAL  CONDUCTIVI  TY  (•.A12.  J 
+  :  ■.2X,  ) 

RE/\13  YKF. 

I  LSEIF  (ANSHK.FO.'N  )  THEN 
GOTO  2^ 

F'l.SE 
GO'TO  2S 
ENDIF 
C 

C  PROVIDE  ANOTHER  REVIEW  OF  EPOXY  ENTRIES 
C 

2"  CALL  CLS 
WRTTE(*.242) 

242  F(2RMaT(;  /.'  you  HAVE  MADE  THE  FOLLOWING  CORRECTIONS  I'O  THE  EP 
+  OXY  EN  TRIES.  ;') 

GO  TO  ,52 
C 

FI.SFIF  (ANSF.EO.  N  )  THEN 
GO  TO  29 
FI.SF 
{>()  TO  .50 
FNDIF 

( 

( 

ckaRACTF.RISTICS********’**’’**** 

c 

c 

CALF  CI  S 
WRTTEC,24.5) 

24.5  F()RMA'T(  '  ■«•»»*■****’■»•*•'  •••••*••»*•»»*♦••♦**•»****»«■»*»••••»» 

^  CHARACTEKISTICS************’ 
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c 

V44  WRi  ri  ('.244)1  1. 

244  FORMAlV  I  N  11  R  C  ARRl!  R  1,1:N(i  I11  (■.A2.  ):  •.2X,  ) 

RI  AD  X('l. 

(' 

C  CARRll  R  1  1  N(j111  .\IL  ST  HI  CiRl'.MTR  lllAN  OR  L:(JUAl,  TO  SUBSTRA  Tt;  LRNG  111 
C' 

IF  (C'l.FI.Sl.j  1(11  N 
CAI.l.  (  LS 
WR1TF(F243) 

PRIN  T  •;  CARRIER  LENGTH  MUS  T  BE  GREATER  THAN  OR  EQUAL  TO  SUBS 
+  I  RATE  LENG  TH. 

GO'T(2  444 
ENDIF 
C 
C 

WRITE(T246)UL 

246  FORMA  T!/, ■  EN  TER  CARRIER  THICKNESS  (’,A2,  ):  •,2X,\) 

READ  FDELC' 

(.' 

write:!  L247)L  K 

24'  I'ORMA  T!  EN  l  ER  C/VITRIER  THERMAL  CONDUCTIVI  TY  (',A12.'):  ■,2X,\) 

READ  *.KC 
C 

.3cS  CALL  (LS 
CW  =  SW 
C 

C  Rl  VltW  CARRll  R  I  N  TRlEiS 
C 

WRTTE(T248) 

248  1  ORMAT(;  ,•  YOU  HAVE  MADE  THE  FOLLOWING  ENTRIES  FOR  THE  CARRl 

+l:r.',;) 

PRIN  T  1.)  LENGTH  '.CL 

PRIN  l'  T'  2.)  WIDTH  '.CW 

PRINT  3.)  THICKNESS  '.DELC 
PRINT  4.)  k  '.KC 

PRINT  * 

PRINT  L  CHANGING  WIDTH  WILL.  ALSO  CHANGE  SUBSTRATE  AND  EPOXY  W 
4-IDTHS.' 

C 

WRITE!*, 249) 

249  FORMAT!/,’  DO  YOU  WISH  TO  MAKE  ANY  CHANGES?  !Y  OR  N)  ■,2X,\) 
READ(*,607)ANSC 

(>()'  I-ORMAT(Al) 

C 

(  MAKE  CORRECTIONS  OR  CHANGES  TO  CARRIER  E/N  TRIES 
C 

II  ( ANS(  .1  (J.  Y  )  THEN 
17  CALI  (  LS 

WRTl  (  (*.M)X) 

GTS  FOR.MAK  ") 

433  PRIM  -.  HIT  CURRLM  I  N  TRY  FOR  LI  NCHH  IS  ',CL,  ,UL 
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WRirH(*,250) 

250  FORMAT(/;  WOULD  YOU  LIKH  TO  CHANGE  LENGTH?  (Y  OR  N)  •,2X, 
+  ’) 

READ(*,(>(W)ANS(T, 

<>09  FORMAT(Al) 

C 

11-  (ANSC'L.LO  Y  )  I  HEiN 
WRlTlr.251)UL 

lORMAIi.  F.NH:R  CARRII  R  l.l-NGTH  (■,A2.  ):  'JX.  ) 

RLAD  TCI 
il  (C1„L1.SL)  mi'N 
CALI.  CLS 
WRHL.(L(i!0| 

010  FORMA !'(,//) 

I'RIN  I  CARRIER  LHNG  I 11  MUS  I  HE  GREAI  ER  LHAN  OR  EQ 

+  UAL  Fo  suhs  i  rai  f;  length.' 

PRIM'  I  111  I'RESEN  I  I  N  FRY  FOR  SUHSIRAIF  I.FNG'l'H  IS 
+  '.SI..'  '.Ul 

PRIM  * 

GO  TO  933 
FNDIF 
C 

EI.Sl  lF'  (ANSCL.EQ.'N')  THEN 
IF  tCL.I  I'.Sl.)  THEN 
CALL  CLS 
WRITE!*, 611) 

611  FORMAT!////) 

PRINT  *,'  CARRIER  LENGl  H  MUST  BE  GREATER  THAN  OR  EQ 
+  UAL  ro  SUBSTRATE  LENGTH.  ’ 

PRINT  *.'  THE  PRESENT  ENTRY  FOR  SUBSTRATE  LENGTH  IS 
+  '.SL.UL 

PRIN'l  ’ 

GO  TO  933 
FNDIF 
tjO  lO  39 
ELSE 
(.o  ro  33 
1  NDIF 

(' 

■59  CALL  CLS 

WRl  1F(  *.612) 

012  FORMAF!  ^) 

PRIM  *.  THE  CURRENl  LN  IRY  FOR  WIUIH  IS  ,CW.'  '.UL 

PRINT  *.'  CHANGING  THIS  LNIRY  WILL  CHANGE:  SUBSTRAFE  AND  ' 

PRIN  I  EPOXY  WIDIHS.  FHEY  ARE  Al  l.  EQUAL.  ' 

Vv  '<1TE(*.252) 

2.>2  FORMAT!/,'  WOULD  YOU  LIKE  TO  CHANGE  WIDTH?  !Y  OR  N)  ',2X,\ 

+ ) 

RFAD!*,613)ANSCW 
6' 3  F'ORMA'I'lAl) 

C 

IF  ( ANSCW.EQ.'Y  )  THEN 
WRITE!  *,253)UL 

2\?  FORMAT!','  ENTER  CARRIER  WIDTH  !'.A2.  );  ',2X,  ) 

re: AD  *.CW 
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i:w  =  c\v 

I  [.Sr;iK  (ANSCW.LiQ.'N  )  ITIIZN 
GO  TO  34 
ll.SF, 

GO  1C)  34 
1  NDll' 

C 

34  C'ALl.  C'LS 
\VRITE(*,614) 

(>14  I-ORMAT(,.//0 

PRINT  THE  CURRENT  ENTRY  FOR  THICKNESS  ’.DELC,’  ',UL 
WRITE!  *,254) 

254  FORMAT!  ,■  WOUFD  YOU  LIKE  TO  CHANGE  THICKNESS?  (Y  OR  N) 

+  2X,0 

REAn(*,615),\NSCT 
()15  T'ORMA'l'(Al) 

C 

IT  (.ASSC'l.l-O.'Y  )  Till  N 
WRIT!  (*.2551)UT 

2551  TORMAT(,'  I  N  H  R  CARRIER  THICKNESS  (•.A2,  );  ■,2X,  ) 

READ  *,DE:I  C 

EI.SEIT  (ANSCT.EO.  N  )  THEN 
GO  TO  35 
ELSE, 

GO'TO  34 
ENDIF 

(’ 

35  CALL  CLS 
WRTTE(*.()16) 

(t  1  EO  R  M  AT(  ////) 

PRIN  T  *,•  THE  CURRENT  EN'TRY  FOR  THERMAL  CONDUCTIVITY  IS  ',KC 
+  ,•  -.UK 

WRTTE(*.256) 

256  FORMAT!/, ■  WOULD  YOU  LIKE  TO  CHANGE  THE  THERMAL  CONDUCTIVIT 

+  Y?(Y0RN)  ■,2X,'.) 

RE;AD!*,617)ANSCK 
6H  EORMATiAl) 

( 

IE  (ANSCK.EQ.  Y  )  THEN 
WRITE!  •,257jUK 

25-  FORMAT!  ENTER  CARRIER  THERMAL  CONDUCTIVI  TY  L.A12. 

+  |:  ■,2X.,) 

READ  *.KC 

ELSEIE  (ANSCK.EO.  N  )  THEN 
GO  TO  40 
ELSE 
GOTO  35 
ENDIE 

(' 

(  .AI  LOW  FOR  ANO  THER  REVIEW  OF  CARRIER  ENTRIES 
C 

40  CALL  CLS 
WRITE!  *.258) 

258  EORMAT!//'/.' 


YOU  HAVE  MADE  THE  FOLLOWING  CORRECTIONS  TO  THE  C 


+  ARRiER  1-:N  I'RH-S.  /) 
(}C)1'()  3(> 


LLSl-IF  (ANSC.I  U.  N')  IHliN 
C}t)  I  () 

(iCIK)  3S 

i;noii- 


+  +  +  P HR ATU re  INF Ui'S •♦♦••♦*****‘''***'*‘****^*********’*'**'*‘ 

c 

CALL  CLS 
\VRITE(*,260) 

+  •** . TEMPERATURE  INPUT . . 


(•  INITIAL  IT.MPERAILRE 
(■ 

WRI  ll  i  \2()1)L  1 

:()1  EORMAIC  I  N  11  R  rUE  INI  TIAL  CHIP  TEMPERATURE  (’.AL  ).  .2X,\) 
lU  AD  M  l 
C 

C  UPPER  AMHITN  1  TEMPERA'TURT- 

c: 

WRITE!  *.2<.2)UT 

262  TORMA  T(  ;  I  N  TER  THE  UPPER  SURFACE  AMBIENT  TEMPERATURE  ('.Al,' 
+  ).  ■,2X..) 

READ  rUPRI 
C 

C  LOWER  AMBIENI  TEMPERATURE 
C 

WRITE!  *.263)UT 

263  TORMA  T{/,'  ENTER  THE  LOWER  SURFACE  AMBIENT  TEMPERATURE  (’,A1.’ 
+  ).  ■,2X,\) 

READ  ‘.LWRT 
C' 

(  RIGH  T  SIDE  AMBIENT  TEMPERATURE 
C 

WRITE(*.264)UT 

264  FORMATf/,'  EN  TER  THE  RIGHT  SURFACE  AMBIENT  TEMPERATURE  C.Al,' 

+  ).  ■.2X.  ) 

READ  M<1 


C  LEFT  SIDE,  AMBIEN  T  TEMPERA!  UKE 
( 

WRTTi  r  .26.S)L'  T 
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2(>5  FORMAT!  , ■  FNTFR  I  HF  l.Fl'l  SURFACF!  AMBIHN  F  rFMPFRA  FURF  (',A1,') 

+  .  ■,2X,',) 
read  M:r 

(' 

C'  I  RON  1  AMBIEN  1  I'i  MFERA'l'URE 
(■ 

\VRri't(T266)L:'l' 

2W>  FOR.M.M  (  .  I  N  l'I  R  TIIF;  FROM'  SURFACE:  AMBIEN  F  l  EMI’FRAFURE  C.Al.’ 

+  j.  •,2\,  ) 

RF,A13  M  'l 
C 

C  REAR  SIDI  .A.MHIl  N  I  1 1  MFI  RA  FL  RE 
C' 

\VRFIl!^2()'’)L' 1  , 

2(.7  FORMA'K  '  1:N1ER  FIIF;  RFAR  SURI  ACE  AMBIEN  F  TEMl’FiRA TURE  (’.Al.  ) 

+  .  ■,2X.  ) 

READ  '.B  F  « 

( 

C 

C  REVIEW  FEMFERAFURE  EN  FRIES 
C 

50  CALL  CLS 
WRITE!  *,268) 

268  FORMA  F!///;,  YOU  HAVE  M.XDE  THE  FOLLOWING  AMBIENT  TEMPERATURE 
+  ENTRIES,',/) 

PRINT  *,’  L)  INITIAL  CHIP  TEMPERA'FURE  MT.’  ’,UT 
PRINT  2,)  UPPER  AMBIENT  TEMPERATURE  ’,UPRT.’  ’,UT 
PRINT*,  3,)  LOWER  AMBIENT  TEMPERATURE  ',LWRT,”,UT 
PRINT  *,  4,)  RIGHT  AMBIENT  TEMPERATURE  ’,RT.'  ’,UT 

PRIN  F  *,•  5  )  LEF'F  AMBIENT  TEMPERATURE  '.LT,'  ^UT 

PRINT  *,'  6,)  FRONT  AMBIENT  TEMPERATURE  '.F'E’  ^U'l 

PRIN'F  *,  7,)  REAR  AMBIENT  TEMPERATURE  '.BT,'  ^UT 

(■ 

C 

WRFFE(*,26y) 

264  f'(3RMAF!,,'  WOULD  YOU  LIRE  FO  MAKE  ANY  CORRECTIONS,’  !Y  OR  N)  ' 

+  .2X,  ) 

REAl3(*,(>!8)ANS'r 
618  FORMA  FiAl) 

C 

C’  MAKE  CORRECTIONS  OR  CHANGES  FO  TEMPERATURE  ENTRIES 
C 

IF  !ANST,EQ.'Y')  then 
42  C/\LL  CLS 

WRITE!*,268) 

PRINT  *,•  1,)  INITIAL  CHIP  TEMPERATURE  ■,FF,-  ^U'F 

PRINT  *,'  2.)  UPPER  AMBIENT  TEMPERATURE  ^UPRT,'  ',UT 

PRINT  3.)  LOWER  AMBIENT  FEMPERAI  URE  ^LWRT.’  ’,07 
PRINT  *,■  4  )  RIGHT  AMBIENT  TEMPERATURE  ',RT,'  ',UT 

PRINT  *,'  5,)  LEFT  AMBIENT  TEMPERATURE  ',LT,'  '.UT 

PRINT  *,'  6.)  FRONT  AMBIENT  TEMPERATURE  ’.FT,’  ’.UT 

PRIN  I  *,  7.)  REAR  AMBIENT  TEMPERATURE  ’.BT,’  ’.UT 

C 

WRITE!  *  2682) 

2682  F'ORMAl!  //,  WHICH  TEMPERATURE  WOULD  YOU  LIKE  TO  CHANGE.’’,/ 
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+.  SI  1. 1(1  /I  Ro  niROUGii  si:Vf;n.  .:x.  , 

Rt-ADc.(>ix)i)i:;i,r 

( 

('  C  ORRIX  I  INI  UAL  1  liMPliRA'l  URK 
11-  (Dhi.r.i-o.  T)  rHi;N 

(  Al  l,  CI  S 
\VR1TF(*,:X)) 

270  1()RMAT(  ) 

PRIN  1  INITIAL  CHIP  I  LMPRLATURE  IS  MT.'  '.UT 
\VRriTC,27I)UT 

271  1  ORMAT(.,‘  ENTER  THE  NEW  VALUE.(’,Ai;)  ',2X.\) 

READ  MT 

(' 

C  ('ORRIX'  r  UPPER  AMBIENT  TEMPERATURE 
El.Sl  Il-  iDEL  r.E(T'2')  THEN 
(  Al  l.  (  I  S 
WRriEC,f<2()i 
(i2()  lORMAli  ) 

PRIN  I  *.  UPPER  AMBIENT  TEMPERATURE  IS  '.UPRT.  .U  P 
WRl  1E(‘,27.7)UT 

2M  1()RMAT(  EN  I  ER  1111  NEW  VAI.UE.('.A1.  )  ■.2X,\) 

READ  Mi  PR  I 

( 

(  (  DRREC  r  l.DWE.R  AMBIEN  1  TEMPERATURE 

El  SMI  (Dl  I  r.E(T  .V)  rilEN 
(  ALL  Cl.S 
WRrrfU*.(>23) 

()23  EORMATf  ) 

PRLN  I  C  LCJWEK  AMBIENT  TEMPERATURh  IS  ■.LVVRT,'  '.UT 
WRI  rE(*,2^MUT 

2^.'  lORMAri;,'  ENTER  THE  NEW  VALUE.('.AL  )  '^Xa) 

READ  M.WRT 
C 

C  .  ORRECT  RIGHT  AMBIENT  TEMPRRAI  URE 
ELSEIE  (DELT.EQ.-4  )  THEN 
CALL  Cl.S 
WRI  rE,(  *.(i26) 

(.2()  l()R.MAr(  ,  ) 

PRINT  V  RIGHT  .AMBIENT  TEMPERA  LURE  IS  ,R  I V  '.UT 
\sRrn(C2'’'iur 

2’’'  iORMAH.  EMI  R  THIM  w \ALUI  .(  Ai,  )  .2X.  ) 

READ  VRl 

( 

(  > ORRECl  l.EE'l  AMBIl  N  T  11  MPERAI  URE 

ELSEIE  (DEI.r.EU  .V)  EHEN 
CA1.L  CLS 
WRITE!  •.(i2V) 

(>2^  EORMAIT/) 

PRIN  I  *.  LEE'I  AMBIE;N1  EEMPl  RA  I  URE  IS  .1.1.  .U  1 
WRl  l  EC.2A|U  r 

2A  EORMAli  .-  ENTER  THE  NEW  VALUE  (•.Al.  )  ■,2X.  ) 

READ  M  E 

( 

(  (  ()RRE(  E  ERGN  I  AMBIENT  EEMPERAEURE 

ELSEIE  (DELT.EO.'f'  )  THEN 
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c:ali,  n.s 
wRrii;(*,632) 

(.22  l{)RMAr(.) 

PRIM  ’.  I  R(JN  r  A.MHII  N  1  I  IMI’I  RAIL  RI  IS  I'l  ,  .L  1 
VVR1'11C.281)L'  I 

2S1  K)RMAI(  ,'  I  N  H-R  nil-  N1  W  VAI.Ul-:.(-.Al.  )  ■.2\.  ) 

Rl  Al)  MT 

(■ 

C  CORRRC  r  Ri-  AR  /VMHlLiN  1'  I  LMl’l-RA  I  L  Rf: 

ELSt'IF  (DELr.hQ.'T  )  IHEN 
GAEL,  GLS 
\VRrrE(*,635) 

635  E'ORMAT(///) 

PRINT  REAR  AMBIENT  TEMPERATURE  IS  '.BT,'  ’.UT 
WRITE(*,284)UT 

284  FORMAT!  ENTER  IIIE  NEWVALUE.C.Al,)  •,2X,.) 

READ  ’.BT 

C 

G  MAKE  NO  GH/\NGES 

EI.SEIF  (DEET.EQ.'O  )  THEN 
4201  GALL  GLS 

WRITE!  *,6320) 

6320  I-ORMAT(  ;V,'  YOU  HAVE  DE.GIDED  TO  MAKE  NO  GORREGTIONS!',/. 

IS  PHIS  GORRLGT.'  !Y  OR  N)  ■,2X,  ) 

READ! *,618)  ANSTl 
IE  (ANSTI.EQ.  Y'i  rilEN 
GOTO  51 

l  l.SElF  !ANS  n.E;0.  N  )  I'HEN 
GO'I'O  42 
ELSE 

GO  TO  4201 
INDIE 
G 

ELSE 
GO  lO  42 
ENDIE 
C' 

C  ALLOW'  FOR  /VNOTHER  REVIEW'  OF  TEMPERATURE  ENTRIES 
G 

48  GALL  GLS 

WRITE!  *.635) 

WRITE!  *,285) 

285  FORMAT!'//,'  WOULD  YOU  LIKE  lO  MAKE  ANY  MORE  GTJRRECl  IONS' 
+  OR  REVIEW  TEMPERATURE  ENTRIES.  lY  OR  N)  .2X,  ) 

READl*,618)ANSTL 
IF  IANSTI..EQ.  Y')  FHEN 
GOTO  42 

f'l.SEIF  (ANSTL.EQ.  N  )  I  HEN 
Go  ro  H 
ELSE 
GOK)  48 
I  NDIF 
G 

I  I.SI  IF  I.ANSI  EO.  N  )  mi  N 
Go  ro  M 
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LLSh 
GO  TO  50 
hNDIF 

( 

( 

( 

(' - 

(■ . ' 

C  "  ■ '  * 

C 

(  PROV1DI-;  A  CGRRI  I.AUGN  hi  IVVl  l  N  nodi-  NUMHI  RRS  AM)  MAIRIX  LOCATION 
C 

M  NLM  =  1 

[)()  (>0  1  =  l.NPL  NWIDI 
DO  M  J  =  I.NWIDI- 
JH{NUM)  =  J 
NL'M  =  NUM  +  1 
(>1  CONTINUL 
(lO  C.'ONIINUL 
( 

NUMA  =  1 
('ll  =  0 
II  =  1 

DO  62  I  =  l.NI’I. 

IH(NUMA)  =  H 
CH  =  CH  +  1 
II'  (C'lI.LO.NWIDL)  TIIT.N 
11=11+1 
('ll  =  0 
LNDII- 

NUMA  =  NUMA  +  1 
(-2  COMINUl 
( 

(  I’ROViDI  Al  ri  RNAllVI  111  Al  INPUI  MI  I'IIODS 

(■ 

(  ALL  (  I  S 
WRiri  CJVO) 

^  j^»p^J|*************#*********^*:^** 


S  v  >  >  V  >  »  •>  -V  •>  -k  -V  ->.  »  V  >  +  >  -- 

‘''^’IILAI  INPU  I 


>11  . 

+  HEAT  INPUT  LO  THE  MICROCIRCUH'  OCCURS  ONLY  ON  HIE  UPPER',/, 
+  SUBSTRATE  SUFACE,  HEAT  INPU  T  IS  ACCOMPLISHED  BY  ONE  OF  ,,/, 

+  THE  FOLLOWING  METHODS;  ,,  , 

+  L)  ENTER  AS  A  TOTAL  HEAT  APPLIED  TO  THE,  CHIP, , , 

+  2.)  ENI  ER  AS  AVERAGE  HEAT  PER  UNIT  AREA,',/, 

+  5.)  ENTER  HEAT  NODE  BY  NODI  ,  , , 

+  4.)  NO  HEAT  INPU  T,'.'/, 

+  PLEASE  SIT.FCT  A  NUMBER  ONE  THROUGH  FOUR,  ',2X,  ) 

RI  AD(*,.'504)SELH 

( 

II  (SLLILFO .1  .OR.SELII  LO  '2'.OR.SFLH,EO.'.5'.()R.SL:LH,F;0,'4')  then 
04  WRI  T!  (•,201  )SFLH 
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:sil  FORMAT! YOU  HAVE  SFLliCTl  D  NUMHIiR  '.Al.’  OF  FOUR  AL  I'F.RNAT 
+  IVES.  IS  THUS  IHE  CORRECT  CHOICE?  (Y  OR  N)  ■,2X.O 

READ(*,304)ANSH 
C 

IF  (/YNSH.EQ.  Y')  THEN 
GOTO  63 

ELSEIF  (ANSH.EO.  N  )  THEN 
GOTO  31 
ELSE 

CALL  CLS 
WR1TE(*,290) 

GO  1 0  64 

63  e:ndif 

C 

ELSE 
GO  I'O  3 1 
ENDIF 
C 

C  DEIERMINE.  UNLIS  KIR  lll.Al  INFL  l 
C 

IF  ((ANSN.EO.  E').AND.(Sl  1  H.EO.  2  ))  IHl-N 
UH  =  'Btu/(hr’in^  2)' 

U/MI  =  Biu/hr 

ELSEIF  ((ANSN.E(T  I  ').ANI),(SELH.EO.  r.OR.SELH.EO.'3'))  THI-:N 
UH  =  '(Btu)  /  (hr) 

ELSEIF  ((ANSN.E0.  S').AND.(SELH,EQ.'2'))  LHEN 
UH  =  ’WATTS/Ccm  2)  ’ 

UAH  =  WATFS  ' 

ELSEIF  ((ANSN.EQ.  S  ).AND.(SELH,EQ.  r.OR.SELH.EQ;3’))  THEN 
UH  =  ■  watts  ■ 

ENDIF 

C 

c 

C  ALLOW  FOR  RE-SELECTION  OF  HEAT  INPUT  METHOD  OR  CONTINUE  WITH 
C  INITIAL  SELECTION 
C' 

C  CHOICE  #1 
C 

CALL  CLS 

11-  (SELH.EQ.'r)  THEN 
63  WRI  FE(’.292)UH 

292  l  'ORMAT(,  ,  YOU  1 1  \\  i-  SELECTED  EO  INPU  I  UFA  I  AS  A  TO  I  AL  HEA  I 

-t-  ,  ,  APPLIED  TO  HIE  SURI  ACE.  . 

ENTER  TOTAL  HEA  I  ,\1’PLIED  lO  I  HE  SURI  A('E  ('.A13,  )  .2X,  ) 

READ  ‘.THEAI 
w.  WRITE!  *.293) 

2W3  FORMAT!/,  IS  THIS  THE  CORRECT  ENTRY  (Y  OR  N)  '.2X,  ) 
READC,304)ATH 

(' 

C  MAKE  HEAT  ENTRY  AND  ALLOW  FOR  CORRECIION 

C 

IF  (ATH.EQ.  Y  ''  THEN 
THPN  =  EHEAT/NPL 
PRINT  • 

PRINT  *,■  TOTAL  Hr;AT  PER  NODE  IS  '.TIIPN,'  '.UH 
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(' 

(  l  ll.l,  HHAT  MAI  KIX  \VH11  |)l  SIKI  I)  V,U-U1-S 

(' 

so  1  =  1, NIDI  I  P 
DO  81  J  =  I.NWIDI- 
lUiAldJ)  =  IIIPN 
SI  ('ONTINLI 

SO  C'ONIINUI 

C 

DLSDIF  (ATll.hO.  N  )  I  HFN 
CN\LL  CLS 
GOTO  65 
ELSE 
CALL  CLS 
\VRITE(*,294) 

294  FORMAT(////) 

PRINT  *;  TO  PAL  MEAT  APPLIED  TO  THE  SURFACE  IS  TREAT,' 

+  ',UH 

GO  I'CD 
ENDIF 
C 

C  CHOICE  #2 
C 

ELSEIF  (SELH.l  0  '2  )  HUN 

67  WRnE(*.295)Ull 

2'^5  FORMATI^  ,  YOU  HAVF  SFl  I  ('I'ED  TO  ENTER  I'HE  AVliRAGE  HI;aT  OV 
+  F:R  I  HF, ./. 

+  UPPER  SUBSTRAl  i:  SURl  AC'E. 

+  EN  TER  THE  DLSIRFiD  HLA  l'  INPU  1  ( '.A!?.').  \2X.\) 

READ  •,/\HL.AT 
C 

C  MAKE  EN  PRY  AND  ALl.OW  FOR  CORRECTION 
C 

68  WRITE!  *,296) 

296  FORMAT!//.'  IS  THIS  PHE  CORRECT  ENTRY.  (Y  OR  N)  .2X,\) 
READ!*.304)ANSHA 
C 

IF  !ANSHA.EQ.'Y')  then 
THPN  =  AiIEAP’SL*SW/NPL 
PRINT  • 

PRINT  *.'  TOTAI.  HEAT  PER  NODE  IS  '.THPN,'  '.UAH 

c: 

('  FILL  HEAT  MATRIX  WITH  DESIRED  VALUES 

C 

DO  82  I  =  l.NDEEP 
D(D  83  J  =  I.NWIDF 
HEAT!  I..I)  =  THPN 
S3  C'ON'I'INUF, 

82  CONTINUE 
( 

FI  SEIF  (ANSIlA.l  O.  N')  FHF.N 
(  ALL  Cl  S 
GO  PCD  67 
FLSI-: 

CALL  CLS 
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WRITEC,2‘)4) 

PRINT  AVERAGE  HEAP  OVER  SUBSTRAI  E  SURI'ACE  IS  -.AHEAT, 

+  '.UH 

GO  TO  (i« 

ENDIF 

(' 

(■  CHOK  E 
(■ 

e;i.se:ii  (Se.i.h  i  (j.  V)  i  hen 

^0  WRIEEG.^^ri 

Z'-r  EORMAK  ,  '»()U  HAVl  SEEECiED  EO  ENIER  HE:a1'  NC.OALLY'.;/, 

+  ENTER  hie:  I'oi  ae  nl  mber  oe  node;s  designated  for  heat  input.  ■ 

+  .2X,  j 

READ  c  rO  ENOD 

c 

C  EHIS  IS  NODE  BY  NODE.  GET  NUMBER  OF  EN  ERIES  THEN  LOOP  UNTIL  ALL 
C  ENTRIES  HAVE  BEEN  MAI3E 
C 

C  FELL  USER  MAXIMUM  ENTRIES  POSSIBLE 
G 

IF  (TOTNOD.GT.NPL)  EHEN 
PRINT  * 

PRINT  MAXIMUM  ENTRY  IS  '.NPL 
CALL  CLS 
GOTO  70 
ENDIF 

(' 

('  M/\KE;  EN  ERIES 
(' 

DO  I  =  l  EOENOD 
NC  =  I 
GALL  CIS 

WREEE'.rL402)Nc:,I(.HNOD 

402  I ORMAI  i  .■  EHIS  IS  NUMBER  I.V  OF  .LG  EN  ERIE:S  ) 

VVRITEi  .248) 

29.S  E'ORMAEl  ;  EN  I  ER  EHE  NODE.  NUMHIiR  FOR  HLAE  INPUT.  ■.2X,  j 
RE/\I3  CNN 

IE  (NN.I  O.O.OR.NN.Cj  E.NPL)  EHEN 
Go  ld  TS 
ENDIF 

\VREEEi’,2VV)UH 

2W  EORMAT(/;  EN  I  ER  HEAT  INPUT  C.Ad,  ). 

READ  *.NHEAT 

HEAT(IH(NN).JH(NN))  =  NHEAT 
CONTINUE 

('  PROVIDE  OPPORTUNITY'  FOR  CORRECTIONS  OR  FURTHER  ENTRIES 

c: 

'3  CAl.L  CLS 

\VRlTE(*,40l)TO'ENOD 

40 1  FORMA  E(  .  YOU  HAVE  MADE  ■.I3.-  NODAL  EN  ERIliS.  ) 

WRITE!  L4(K)) 

4(H)  FORMAEI  .  DCJ  YOU  WISH  EO  MAKE  ANY  MORE  EN  ERIES  OR  CORRECTIO 
-t-NS.’  (Y  OR  N )  'JX.  ) 


IF  (AJIN.EQ.'Y  )  THEN 
CAI.L  CFS 
GO  rO  70 

EEsi:iF  (ahn.eu.  N  )  then 
GO  TO  72 
El.SE 
GOTO  "7 
72  ENDIF 
G 
G 

G  NO  HHA  I  INFU  l  -  GIIOIGH  #4  -  HEAT  MATRIX  STAYS  INITIALIZED  AT  ZERO 

G 

LI.SEII-  (SIT.II.EO.  4  )  niEN 
GO  TO 
I  NDIF 
G 

«()  GONTINUL 

G 

G 

G 

G  DETERMINE  INCREMENTAL  MFASUREMEN  LS  IN  HIE  X  AND  Y  DIRECLIONS 
(' 

DELX  =  SLNDEEF 
DEFY  =  SWVNWIDE 
G 

C****-'-* . GENERATE  CONST/\NTS  FOR  THE  SUBSTRATE  LAYER . 

G 

C  LER-  OR  RIGHT  EDGE  TO  OUTSIDE 
SYLR  =  2  *  KS  *  DELX  ’  DELS  /  DELY 
G 

G  FRONT  OR  BA(  K  EO  OU  TSIDE 

SXFB  =  2  •  KS  ’  DELY  *  DELS  DELX 
G 

G  INNER  MATRIX  MOVEMENT  IN  THE  Y  DIRECTION 
SYY  =  KS  '  DI  LX  '  DELS  DELY 
G 

G  INNER  MATRIX  MOVEMENT  IN  THE;  X  DIREGTION 
SXX  =  KS  '  DELY  '  DELS  DELX 
G 

(  SUBSTRATE  TO  TOF  SURTAGE 

SZ  T  =  2  *  KS  *  DELX  ‘DELY  DELS 

G 


(  SUBSTRATE  T(J  EFOXY 

SZE  =  2  •  DELX  •  DELY  ({DELS  KS)  +  (DELE:,KE)) 

G 

CONSTANTS  FOR  EPOXY  i.aYER****’********* . 

G 

G  LEF'T  OR  RIGH  T  EDGE  TO  OUTSIDE 
EYLR  =  2*KE*DELX*DELF/DELY 
G 

G  FRON  T  OR  BAC  K  EDGE  TO  OU  TSIDE 
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HXFH  =  2*Kh'DI  l,Y*DFLlLDia,X 


INNI  R  MAIRIX  NKA  l-MFN  l  IN  IHl  DIRI  (’HON 
IVY  =  Ki'Di  1  x*di;lf  DHLY 

l  l'OX'i  lOCARRlFR 

i  /(  =  l.X’Ol  l.Y  ((01L1-/R1'.)+(DL-1.(:  KC)) 

INNFR  MAI  RIX  MOVEMFN  I'  IN  TIIH  X  DIRECTION 

Exx  =  ke*di.:ee*dely,delx 

•••'•’’’‘’■'■’“’■’‘’GENERATE  CONSTANTS  FOR  INTERIOR  CAJ^RIER  LAYER* 

I, LIT  OR  RIGHT  EDGE  TO  OUTSIDE 
CYLR  =  2*KC*DELX*DELC.DELY 

INNER  MATRIX  MOVEMENT  IN  THE  Y  DIRECTION 
CYY  ^  KC’DEI.X’DELCDEEY 

INNER  MATRIX  MOVEMENT  IN  THE:  X  DIRIT'TION 
CXX  =  K(  ’D1  EC’DELY  DELX 

CARRIER  lO  liOITOM  OU  TER  EDGE 
C/Ii  =  :’K(  Dl  l  X’DE.EY  DEE( 

(  ARRIl  R  EO  ER(JN  I  EDGE  It  NO  EAR  EXISTS 
CXE'B  =  :*KC’DEI,Y’DEECDELX 


■•'COEEEK  U  N  IS  TOR  SUBSTRATE.  EEOXY,  AND  INTERIOR  CUVRRIER  LAYERS* 


DO  90  I  =  l.NTI. 
N  =  1 

IB  =  NFE  +  I 
ID  =  2’Nin +1 


IT  ((IH(I).E0.1.OR.IH(I).EQ.NDEEP).AND.(JH(n.EQ.l.OR.JH(T).E(^. 
+  N\VIDE))  THEN 

DE  TERMINI  CONNECTIONS  EOR  TOP  PAYER 
IE  (HEAT(IH(I).JH{I|).E;CTO.O)  Till  N 
NC’ONH.N)  =  (1 
EESE. 

NCON(EN)  =  ' 

I  NDIT 


CONNIr  TIONA  TOR  EPOXY  LAYER 
N(  ONdB.N)  =  h 


76 


c 

(■ 

('  Lt  l'l  AND  RIOH  1  C'Ol-MK'll  N  IS  Dl  IM  NDINCi  ON  WHICH  l-.i:)C;i: 

(  lil'lllX.I 

II  (JHiH.l  U.l)  HlI  N 
(■  1  I  I'l'  COM'I-ICII-N  1 

COI  ImI.N)  =  LK 
CtM  I  HH.N)  =  I  'll  K 
N  =  N  +  1 
NC(JN(l.Ni  = 

N('()N(IB.N)  =  -.r-Si 
C  K1C;H  1  COl'M'lCIIN  1 

col  l-il.N)  =  SVA 
COl-l'(IB.N)  =  I'VY 
N  =  N  +  1 

NCONd.N)  =  U)’(I  +  1)+1 
NCON(IB.N)  =  10*{IB+1)+1 
(  RIGHT  TDGR 

Id.SIill'  (JIUD.HQ.NWIDli)  THIiN 
C  I.[:FT  COI  FFICir-NT 

COHFd.N)  =  SYY 
COiddlB.N)  =  FYY 
N  =  N  +  1 

N(’ON(FN)  =  10’(I-U  +  1 
NC(JNdn.N)  =  lOdlB-H  +  l 
C  RK)H1  COFFFlCUiN  r 

COld(I.N)  =  SYFR 
COl-.FdB.N)  =  FYl.R 
N  =  N  -  1 
N{  ON(  I.N  )  =  'Ml 
NCONdl'.N,  =  'S41 
INDII- 

( 

C  I  RON!  AM)  BAC  K  coin  II  II  NFS  1)1  I’l  NDINO  ON  WHICH  1:1)01 
(  I  RON!  FDCH- 

II'  dlld), 1:0.1)  niFN 
(  I  RON  I  ( OFinCll  N "1 

col  |■■d,N)  =  SXI  B 
col  l  (IB.Nj  =  F.Xl'B 
\  =  N  +  1 
NCONil.N)  =  7.s:i 
NCON'dB.Nj  =  7.^21 
(  BACK  COFFFICIFNT 

COFl-(I.N)  =  SXX 
COEFdB.N)  =  FXX 
N  =  N  +  1 

NCONfl.Nj  =  lOdl  +  NWIDld+l 
NCONdB.N)  =  10*(IB  +  NWII)1:)+1 
C  BACK  F'DOF 

FI  SI  (F  (IHdj.FO.NDF.hld  IHFN 
(  IRON  1  COFl'l'IClFlNI 

(-OI:l'([.Ni  =  SXX 
( OFl'dB.N)  =  FXX 
N  —  N  1 

NCONd.N)  =  Dt'd-NWIIJF  !+ I 
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NCONdH.N)  =  10dlB-NVVlDH)  +  l 
C  BACK.  COHFFICIl-N  I 
COfd'(I.N)  =  SXI'B 
C'OFEdIB.N)  =  FXr-B 
N  =  N  +  1 
NX'ON(l,N)  =  75.^1 
NCONdH.N)  =  ".XM 
ENDII- 
C 

C  r01>  COI  J  I  ICII-N  T 

COl  F(I.N)  =  SZT 
COIl(IB.N)  =  S/l- 
N  =  N  -  1 
NCONd.N)  =  '.Ml 
NCONdli.N)  =  l(ld+  1 
(  BOl'IOM  COFld-lClFN  I 
COl  I'd. N)  ^  S/F. 

(  (d  l  ilB.N)  =  F/( 

N  =  N  +  1 

N(  ONll.Ni  =  lOdl  +  M’l.l+l 
NCON'dB.N)  =  10*dB  +  Nl>l,)+l 
(  Ill-.AI  INPL  1 

II  dii:Ai(ind),Jiid)).NF..o.O)  iiifn 
COId-d.N)  =  HF.ATdHdl.JHd)) 

N  =  N  +  1 
N('ONd.N)  =  VVVI 

I  N  I)  1 1 
C 

C 

( . RONT  AND  BACK  FIDGES  EXCl.UDING  CORNERS*’'*****’ 

C 

C 

M  SI  II  (dII(I).I-Q.1.0R.IHd).EU.NPEF,P).AND.(,md).NF.1.0R.Jlld) 

+  .NI  NWIDI  ,,  rUHN 
C 

(  Dl n  RMINI-  NL'MBIIR  OF  CONNECTIONS  FOR  SUBSTRATE;  I.AYI  R 

II  dii;Ardiid).Jiidi).i  Q,().(i)  iiifn 
NC  ( )N(  l.N )  =  () 

MM 

NCONd.N  (  =  7 
1  NDIl 

C 

(  Dl  ri  RMINE  NUMH1:R  of  CONNFCnONS  FOR  IT’OXY  I.aYFR 
NCCTNdH.N)  =  0 

( 

(  I.FFT  COEFFICIFNl 

( OF  Fd.N)  =  SYY 
COEF{IB.N)  =  EYY 
N  =  N  +  I 

NCONd.N)  =  l()*d-l)+l 
NCONdH.N)  =  l()*dH-l)  +  l 
(  RIGHT  COEFFICIENT 
( OEFd.N)  =  SYY 
COFFdH.N)  =  EYY 
N  =  N  +  1 
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NCONil.Ni  =  1 0*  1 1  +  1 1 -I- ! 

NCONtlB.N)  =  l()Mlli+i)+l 

( 

C  1  ROM  ANO  HA(  k  t  ol  l  I  It  II  N  IS  DI  FI  NDI  N  T  ON  WHICH  i  nOl- 
t  1-RON  r  I-.DCjI- 

11-  ,!H(i).i-:u.i]  riiRN 
C  I  RON  T  C0R1'1'1C11:N1' 

COl  l-(I.N)  =  SXFB 
COFF(IB.N)  =  FXFB 
N  =  N  +  1 
NCON(I,N)  =  7521 
NCQN(IB,N)  =  7521 
(  BACK  COFFFICIFNT 

COF.F'(FN)  =  SXX 
cof;f(1B.N)  =  Fxx 
N  :=  N  +  1 

N(  ()N(FN)  =  10*(1  +  NW1HF)+1 
NCON(IB.N)  =  10’(IB  +  NWIDF;)+1 
C  BACK  l  OtlF 

FI.SFII-  (111(1). FO.NDI  FF)  IHFN 
(  I  RON  1  coi  l  I  ICl!  N  I 
COFF(I.N)  =  SXX 
COl  FdB.Nj  =  FXX 
N  =  N  +  1 

NCONH.N)  =  1()iI-NWI13F:;  + I 
NCONHB.N)  =  H)'(1B-NA  11)1  )+  I 
(  BACK  COl  1  l  ie  II  N1 

( OFFd.N)  =  S.XFB 
COF.FdB.N)  =  F.XFB 
N  =  N  +  1 
NCON(l,Nj  =  75,d 
NCONdB.N)  =  7531 
FNDIF' 

( 

C  FOR  COld  'FlClFN  r 

coFFd.N)  =  szr 
coiFdB.N)  =  szf; 

N  =  N  +  1 
NCONH.N)  =  7511 
NCONHB.N)  =  lOd+1 

( 

C  BoriOM  COFFFICIFN  1 
COFFd.N)  =  SZF 
COFFHB.N)  =  FZC 
N  =  N  +  1 

NCONH.N)  =  Unl  +  NI’l.j+l 
NCONHB.N)  =  li)-HB  +  NFl.)+ 1 

( 

(  litAllNFLI 

IF  HlFAldllH,.. 111(1)). NF.O.O)  IIIFN 
(OFFd.N)  =  HFAldlliF.JHHi) 

N  =  N  +  I 
NCONH.N)  =  WI 
FNDIF 

( 
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(' 

(''‘*”*******LEFr  /VND  RIGHT  EDGES  I  XCL.UDING  CORNERS******^*’*******’***** 

G 

C 

I  LSI- It-  ((.III(Ii.I-(J.1.C)R..I11|I).I  O.N\MDI-|.AND.(III(1}.NI-.1.GR,1H(1) 

-r.M-.NDEEPil  THI  N 

(■ 

(  Dl- n-RMINI  NL  MRl  R  Ol-  ( ONNiiCl  IONS  1  OR  SUliS  I  RA  lT  L.AYIIR 
II-  (llEAI(lll(li„)ll(l)).Nl  ,0.0)  IIIEN 
NCON(l.N)  =  ^ 

ELSE 

NCON(I.N)  =  (1 
ENDIF 

(' 

C  DE  EERMINE  EHE  NUMBER  OF  CONNECTIONS  FOR  EPOXY  AND  CARRIER  LAYERS 
NCON(IB.N)  =  ft 

NCONfID.N)  =  ft  ( 

C 
C 

C  L.EI-T  AND  RIGHT  COEI  EICIEN  IS  DEPENDING  ON  WHICH  EDGE 
C  LEET  EDGE 

IF  (JHfl).EQ.l)  THEN 
C  LEFT  COEFE'ICIENT 

COEF(LN)  =  SYLR 
COEF(IB.N)  =  EYI.R 
COEFdD.N)  =  CYLR 
N  =  N  +  1 
NCON(I.N)  = 

NCON(IILN)  = 

NCONdD.N)  = 

C  RIGH  T  COEl-FICIEN  T 

COEI-'d.Nj  =  SY^ 

COEE'dli.N)  =  EYY 
COEFdD.N)  =  (  YY 
N  =  N  +  1 

NCONd.N)  =  10*fl-t-l)-t-l 
NCONdli.N)  =  10*dB+l)+l 
NCONdD.N)  =  10*(1D+1)  +  1 
C  RIGHT  EDGE 

ELSEIF  (JHd).EQ.NWIDE)  THEN 
C  LETT  COEFFICIENT 

COEF(LN)  =  SYY  ' 

COEFIIB.N)  =  EYY 
COEF(ID.N)  =  CYY 

N  =  N  +  1  '■ 

NCON(LN)  =  10*(I-1)+1 
NCON(IB,N)  =  lO*dB-I)+l 
NCONdD.N)  =  I0*(ID-1)  +  1 
(  RIGHT  COEFFICIENT 

COE:-'d.N)  =  SYLR 
COIFdlLN)  =  lYLR 
( OE.I-dD.Nj  =  (  YLR 
N  =  N  +  1 
NCONd.N)  =  'Ol 
NCONdB.N)  =  'Ol 
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U  o 


NCON(ID,N)  =  7541 
ENDIF 

(' 

C  FRONT  COEFFIC:iFNT 
COEF(l.N)  =  SXX 
C:OEF'(IB.N)  =  EXX 
C'OEFlID.N)  =  CXX 
N  =  N  4-  1 

NCONd.N)  =  10‘(I-NWIDE)  +  1 
NCONdITN)  =  1()*(1B-NWIDE)+1 
NCONdD.N)  ==  10*(in-NWIDI')+l 

(' 

C  BACK  C  OEI  FICIFN  1 
t:OEF(I,N)  =  SXX 
COEFdB.N)  =  FXX 
C'OEFdD.N)  =  CXX 
N  =  N  +  1 

NCON(I.N)  =  10*(I  +  NWI01)+1 
NCON(IB,N)  =  10*(IB  +  N\V1DE)+1 
NCON(in.N)  =  1()*(10  +  N\V10E)+1 

l  OP  CQEFI  ICIEN  I 
C0EF(I,N)  =  SZT 
CQEF(IB,N)  =  SZE 
COEF'dD.N)  =  EZC 
N  =  N  +  1 
NCON(I,N)  =  7511 
NCON{IB,N)  =  10»I+1 
NCON(ID,N)  =  10*IB+1 
C 

C  BOrrOM  COEFFICIENT 
COEF(I,N)  =  SZE 
COEF(IB.N)  =  EZC 
COEFdO.N)  =  CZB 
N  =  N  +  I 

N('ON(I.N)  =  10*(I  +  NPI,)+1 
NCONdB.N)  =  lOdlB  +  NPI  )+l 
NCONdD.N)  =  '5(.l 
C 

(  IIIATINPLI 

IF  dIEArdlid),Jlld)).NT.(M))  IdFN 
COEI(I.N)  =  IIFATdlId).,llld)) 

N  =  N  +  1 
NCONd.N)  =  9SW1 
ENDIF 
C 
C 

C  ’’’•••DETERMINE  COEFTICIENTS  FOR  ALE  NODES  NOT  TOUCHING  /\N  EDGE . 

C 

c: 

EI.SEIF  (dHd).NE.1.0R.IHd).NE.NDEEP).AND.(JH(I).NE.1.0R.JH(I) 

+  .NE.NWIDE))  THEN 
C 

C  DETERMINE  NUMBER  OF  CONNECTIONS  FOR  SUBS  I  RAI  '-;  LAYER 
IF  diEArdHd),JHd)).NE.n.n)  niEN 
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NCONd.N)  =  ■ 
tiLSl- 

NC:()N(l,Nj  =  () 

HNOIF 

(' 

(  DFIHRMINI-;  NUMBl-R  Ol-  CONNliCHONS  FOR  liPOXY  AND  CARRiDR  l.AYDRS 
NCONiIB.N)  =  6 
N('()N(ID,N)  =  b 

(' 

C  LEFT  COEFFICIENT 
COEF{I.N)  =  SYY 
COEF(IB.N)  =  EYY 
COHF(ID.N)  =  CYY 
N  =  N  +  1 

NCONd.N)  =  lO’d-d  +  l 
NCON(IM.N)  =  10*(IB-n+l 
NCON(ID,N)  =  10*(ID-1)+1 
C 

C  RIOIFF  COl  FF'ICIEN  F 
COI  F(I,N)  =  SYY 
COFF'dB.Nl  =  I  Y\ 

COEFiID.N)  =  CYY 
N  =  N  +  1 

NCONd.N)  =  l()*(i  +  lj+l 
N(’ONdB.N)  =  11)*(1B+Il  +  1 
NCONdD.N)  =  1()*(ID+1)+1 
C 

C  FRON  r  COEFFICIEN  r 
COEF'd.N)  =  SXX 
COEF(IB.N)  =  i:XX 
COEFdD.N)  =  CXX 
N  =  N  +  1 

NCONd.N)  =  10*(I-NWIDE)  +  1 
NCONdU.N)  =  10*dB-NWlDE)+l 
NCONdD.N)  =  10*(ID-NWIDE)  +  1 
C 

C  BACK  COEFFICIENT 
COEF'd.N)  =  SXX 
COEIdIB.N)  =  EXX 
COEFdD.N)  =  CXX 
N  =  N  +  1 

NCONd.N)  =  10*d  +  NWIDE)+i 
NCONdB.N)  =  1U*(1B  +  N\V1DE)  +  1 
NCONdD.N)  =  11)*{ID  +  N\VIDE)+1 

( 

{  FOP  CCJEFFICII  NI 
COEFd.N)  =  S/d 
COEFflB.N)  =  S/F 
coEFdD.N)  =  f:/c 
N  =  N  +  1 
NCONd.N)  =  7S1 1 
NCONdB.N)  =  10*1  +  1 
NCONdD.N)  =  10*IB+1 

( 

(  BorrOM  COFTMCIEN'r 
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COtF'd.N)  =  S/I- 
COHFdH.Ni  =  \  VC 
COi-FdD.N)  =  C/li 
N  =  N  +  1 

NCONd.N)  =  ll)'d  +  M'l.)-r  1 
Nt'ONdllN)  =  Ki'ilH  +  SPI  i-t-l 
NCONdllNi  =  'Sdl 
C' 

(  Ill-AllNPLl 

II-  dii-.-\rdiid)..iiid)).Nd-.o.())  1  iii-N 
C'oi  1  (I  N)  =  111  Aidiidj.Jiid)) 

N  =  N  +  1 
NCONd.N)  =  W91 

i;ndif 

c 

FNDIF 

90  CONTINUF 

C' 

C 

C****‘***COFFFlCIFNTS  FOR  CARRIFR  FAYFR  FRONT  AND  BACK  EDGES . 

(' 

(  DF  ri-RMINF  ('ONSTAN'IS  FOR  CARRIER  COEFFICIENTS 

C 

(  FAR  SI/F 

I  AR  =  (CF-SI.)  : 

C 

("*"”*•••**'■•**1  AR  coi  l  FK  ll-N  IS  IF  NECCIfSARY**"*""-**********”-**'**** 

^  ■  «  V  ^  «  n  .n  «  «  ^  V  C  ^  ^  *  i<  ,K  *  *  (;  ■».  x  ■  s  .  v  t  I*  ^  ^  -v  V  »  i,  *  .<1  +  »  »  * 

c 

IF  iFAR.NI-.O.O)  nil  N 
C 

(  DF  IFR.MINF  III!  N I  C  C  I  SAR Y  ('ONS 1  ANl'S 
C  CONSI  AN  IS  FOR  EAR  NODES 
I  I  I  FT  OR  RICH  I  l  O  FX  ri  RIOR 

CYFRF  =  :*kc*f;ar*defcdffy 

C 

(  I  1  F’l  OR  RIGIM  FO  IN  FFRIOR 

c  YYF.  =  I  ar*kc*df.fc.df:fy 

( 

C  I  RON  !  OR  BACK  FO  EXTERIOR 
CXFBF  =  2*DFFY*KC*DEFCEAR 

( 

(  FRONT  OR  BACK  TO  IN  FFRIOR 

CXXF  =  DFd,Y'DFFC*KC'2d:,AR  +  Dl  EX) 

( 

(  lOPORBOEFOM 

(  /I  AR  =  K(  •EAR*F)I  EY’2  DFE( 

( 

(  Dl  I  ERMINE  I  AR  (  OI  I  I  RUN  IS 
DO  1(H)  I  =  I.2'NWIDI 
II  =  VNPI  +1 
N  =  i 
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('  l-.AR  NODI-;  CONNFCTIONS 
NC'QN(Il-,N)  =  6 

(' 

('  (ORNFRS 

IF  ((IH(l).F;0.1.OR.IH(I).EQ.:).AND.fjn(!).F;0.1.()R.JlI(r).F0. 
+  N\VIIDF.))  TnF;N 
('  Fl  l'l  I  IXiF 

II-  (JH(1).F;u.1)  IHFN 
(  FFl-'l  C'OFl  l-ICIFN  r 

COl  1  (IF.N)  =  CYI.RF 
N  =  N  +  1 
N('()N(11-,N)  =  :’5,S1 
(  RKill  1  C  ()F;1FIC'11-.N  I 

('C)I  l■'(IF,N)  =  CYVF 
N  =  N  +  1 

NCONdl.N)  =  1()*(IF+Ij+1 
Fl.Sldl-'  (JH(l).FaNWlDF.)  TllFN 
(  RKiHlFIXiF 
('  1,1  1-T  COl  I  I  IC'lFN  l 

cof:f(1i;,N)  =  c'yye 

N  =  N  +  1 

NCON(IE,N)  =  10*(IE-1)  +  1 
C  RIOin  COEFFICIENT 

(.:OFF(IE,N)  =  CYLRE 
N  =  N  +  1 
NCON(IE,N)  =  7541 
FNDII' 

C 

C  I  RON  I ,  HACK  AND  TOP  COEFFICIENTS 
('  FRONT  EDGE 

IF  (IIIdl.EO.!)  rHF;N 
C  F'RONI  COEFFICIENT 

COl  F'dE.N)  =  CXFHE 
N  =  N  +  1 
NCON'dF.N)  =  1 

(  !iA(  K  COEFFICII  N  I 

COEF(IE.N)  =  CXXE 
N  =  N  +  1 

NCONdl  ,N)  =  10*(:*NI'I  +.IIId))+l 
(  lOP  (  ()l  1  1  ICII  N  I 

COl  1  (IF.N)  =  C/F.AR 
N  =  N  +  1 
NCON(IF.iN)  =  7521 
(  HACK  EDGE. 

FLSEIF  (111(1). E0.2)  THEN 
C  FR(7NT  COEFFICIENT 

COEFdE.N)  =  CXXE 
N  =  N  +  1 

NCON(IE.N)  =  U)*(.1*NPL-NWIDE+JH(I))+1 
C  HACK  COEFTICIENT 

COEFdE.N)  =  CXETSE 
N  =  N  +  1 
NCON'dF.N)  =  7.S31 
C  FOP  COEFTICIENT 

COEF(IF-.N)  =  C/EaR 
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N  =  N  +  1 
NC'ONdF.N)  =  7531 
INDIl- 

(■ 

(  BO  i  rOM  COt-.FFlClLNT 

CQFFdlF.N)  =  C/.F/XR 
N  =  N  +  1 
NCONdL.N)  =  ^561 
C' 

C  l  AR  NC)[5FS  ON  FRON  T  FDGF  FXCLUDING  CORNFRS 

liLSFlI-  ((IH(I).EQ.l).AND.(JH(l).NFd.OR.JII(l).NE.NWIDE))THFN 
C  I.EFT  COEFFICIEN'T 

COEFdE.N)  =  CYYE 
N  =  N  +  1 

NCON(IE,N)  =  l()*(3’^NT!.+Jd(I)-l)+l 
(■  RKJITTEDGE 

(■0FF(H:,N)  =  CYYF 
N  =  N  +  1 

NCON’di:,Nj  =  11)'{3’'NPI  +,llld)  +  l)+l 
(  TRON'T  C Old  FK  II  N  1 

('OFl(IF..N)  =  (XTUT 
N  =  N  +  1 
NCON'dli.N)  =  7521 
(  BA(  K  COFFFICIFN  T 

COlFdE.N)  =  CXXl 
N  =  N  +  1 

NCONdE.N)  =  lU*(2'NPL+JHd))  +  l 
('  TOP  COEF  FICIENT 

COEFdE.N)  =  CZEAR 
N  =  N  +  1 
NCON(FE.N)  =  7521 
C  BOri'OM  COEFFICIENT 

COEF(IE.N)  =  CZEAR 
N  =  N  +  1 
NC;ON(lE,N)  =  7501 
C 

(  F  AR  NODES  ON  HACK  EDGE  EXCLUDING  CORNERS 

Id,SFIF((IfId).E0.2).A.ND.(JH(I).NEd.0R.JHd).NF.NWIDE))THEN 
C  I  F  I  T  COFFl-ICIENT 

COEFdE.N)  =  CYYE 
N  =  N  +  1 

NCONdE.N)  =  IU'(3*NP1,  +  NW1DF:+Jlld)-I)  +  I 
C  KICilll  COFFFICII  N  T 

COEFdE.N)  =  (  YYF 
N  =  N  +  I 

NC'ONdl-.N)  =  10*(3'‘NP1  +N\VID1  +IFFd)+l)+l 

(  FRONT  cof;fficifn  1 

COEF(IE.N)  =  CXXF- 
N  =  N  +  1 

NCONdF-.N)  =  U)N3*NPI,-N\VIDF+JIFd))+l 

(  Hack  coff  ftcien  i 

( OIFdE.N)  =  CXF-dU 
N  =  N  -»■  1 
NCONdE.N)  =  7531 
(  TOP  C'OFIT-ICIFNT 
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CObFdE.N)  =  C/FAR 
N  =  N  +  1 
NC:i)N(IF.N)  =  "Ad 
(  HO  FIOM  COFI-l  R  UM 
(OFl'dF.Nl  =  C'/FAR 
N  =  N  +  1 
NC'ONdF.N)  =  "Mil 
FNOn- 

100  CON'llNL'F 
( 

C 

C 

,,|TRQN1'  and  BACK  FDGF  COEFFICIENTS  IF  AN  EAR  EXISTS*** . 

c: 

C 

c: 

DO  101  I  =  l.NPL 
N  =  1 

IC  =  2*NPF  +  I 


C 

C 

C* 

C 


V  ¥  *  :*  •*; 


‘CORNERS'''* 


II'  (dHd).EQ.l.OR.IHd).F0.NDEFP)  AND.(JH(I).E0.1.ORJIld) 
+.f,(tn\vidf;))  riiFN 
c 

C  NUMBER  CM-  CONNI  (' 1  IONS  PI  R  NODE, 

NC{)N(IC.N)  =  0 
C 

(  l.I  Id  AND  RICH  I  COE.I  I  IC  ILN  IS  Dl  PENDING  ON  WHICH  EDGE 
C  I, Eld  EDGE 

II-  (JHd).EQ.I)  rilEN 

c  i.EFi  coei  ficie;nt 

CCJEl-dC.N)  =  CYER 
N  =  N  -f  I 
NCON(IC,N)  =  7,AM 
C:  RIGHT  COEmCIENT 

COEF(IC.N)  =  CYY 
N  =  N  +  1 

NCON(IC.N)  =  10‘(IC+1)  +  1 
(  RIGHIEDGE 

I  ESEIF  (JH(I).ECTNWIDE)  I  HEN 
C  I, FIT  COEFFICIENT 

COEF(IC,N)  =  CYY 
N  =  N  +  1 

NCONdC.N)  =  10*dC-l)-t-l 

c  RIGHT  coefficie:nt 

COEF(IC.N)  =  CYER 
N  =  N  +  1 
NCONdC.N)  =  ">41 
1  NDIl 

( 

C  I  RON  I  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  E.DGE 
(  IRON r EDGE 


II  (IIIdi.ECJ.l)  nil  N 


< 
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C  FRONT  COFFFICIENT 

COEF(IC.N)  =  CXXE 
N  =  N  +  1 

NCONdC.N)  =  10‘(3‘NFL+JH(I))+1 
C  HAC  K  C'OIT  I  K'lEN'l 

C'OFFdC.N)  =  C'XX 
N  =  N  +  1 

NC:ONdC',N)  =  10dIC'  +  N\VlDE)+l 
(  HAC  K  IDCiF 

11  SI  !1  dlldi  l  (TMJI.I  Id  nil-.N 
(  I  KOM  C  C)1  I  1  1C  IIM 

C'OI-IhIC'.N)  =  C'XX 
N  =  N  +  1 

NCONdC.N)  =  lOdlC-iNVVlDFj  +  l 
C  HACK  COFl  1  K:I1  NT 

COEFdC.N)  =  CXXE 
N  =  N  +  1 

NCONdC'.N)  =  lOd.^^NPL  +  NWlDE+JIldll+l 
I  NDll 
C 

c:  lop  coefficif:nts 

C:OEF(IC,N)  =  EZC 
N  =  N  +  1 

NCON(IC.N)  =  U)*(IC-NPL)+1 
C  HO  ITOM  COEFFICIENT 

COF.I'(IC,N)  =  CZB 
N  =  N  +  1 
NCON(IC,N)  =  7.S61 
C 
C 

‘ . *****FRON  r  AND  HACK  EDGES  EXCLUDING  CORNERS*****'*****’ 

C 

C 

11  SI  IF  (dlld).E0.1.OR.IIId).EQ.NDEE;P).AND.(JIId).NE.l.{JR. 

+  111(1). NF.NWIDI  ))  HI!  N 
C 

C  NCMHER  OP  C'ONNl-CnONS  PER  NODE 

NCONdC.N)  =  (. 

C 

C  1,1  PT  COI-PPICIP.N  1 

CC^EFdC.N)  =  CYY 
N  =  N  +  1 

NCONdC.N)  =  10*dCd)+l 

c 

C'  RIGHT  COEFFICIENT 

COEFdC.N)  =  CYY 
N  =  N  +  1 

NCONdC.N)  =  10*dC+l)  +  l 
C 

c:  I  RON  1  AND  HACK  COEPTICIEN  PS  DEPENDING  ON  WHICH  EDGE 

C  FRONT  EDGE 

IF  dH(I).EO.l)  THEN 
C  FRCJNT  COEPTICIEN  I 

C'C.)EFdC:.N)  =  CXFBE 
N  =  N  +  1 
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NCONdC.'.N)  =  li)d3*Nl'l.+.Ill(l)j+l 
(  BAC  K  c:oi-,i-Fi('n:N  1 

C.’OFI'dC.N)  =  CXX 
N  =  N  +  1 

NC:C)N(IC:.N)  =  10*(IC+NWIDF)+1 
C  BACK  FDGH 

l-l.Sf-IF  (IH(i).HO.NDHF.P)  THF.N 
C'  I  RON'l'  COI  I'I  ICIFNT 

COII(IC.N)  =  CXX 
N  =  N  +  1 

NCCN(IC,N)  =  10*(1C-NW1DL')  +  1 
('  BACK  COFFF'ICll  N  1 

coFFdc.N)  =  cxfbf: 

N  =  N  +  1 

NC'CJNdC.N)  =  l0d3’Nl>L  +  NWIDF:+JH(l))+l 
FNddl- 

(  ( 

(  roi'  ( (H  id  icn-N  I 

C(HddC.N)  =  1/C 
N  =  N  +  1 

NC'CNdC.N)  =  l()dlC-NI’F)+l 

c 

C  BOFIOM  CCH  1  1  K  IFN'l 
COF.l'dC.N)  =  C/B 
N  =  N  +  1 
NCONdC.N)  =  7561 
F.N'DIl 

101  CONTINUE 

C 

C’****l  RONT  AND  BACK  EDGE  CARRIER  COEFFICIENT  IF  NO  EAR****** . 

C 

ETSEIF  (EAR.I'(TO.O)  THEN 
DO  102  I  =  l.NPL 
N  =  1 

ic  =  :'NPi.-t-i 

( 

c 

II  (dlld).F.O  1.0R.IIIdi.FU.NDEF:P)./\ND.(JIId).F:(T1.0R..Illd). 

-^FO.NWIDI  ))  nil  N  i 

( 

(  NUMBER  OF  ('ONNFCTIONS  PER  NODE 

NCONdC.N)  =  6  * 

( 

(  I.EIT  AND  RIGHT  COEFEICIEN  I  S  EH  PENDlNC.  ON  WHICH  EDGE 
(  I. EFT  EDGE 

11  (JH(Ii.EO.l)  THEN 
(  I  EFT  COEFFICIENT 

COl  IdC.N)  =  CYLR 
N  =  N  +  1 
NCONdC.N)  =  7551 
C  RIGHT  COEFFICIENT 

COFF'dC.N)  =  CYY 
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N  =  N  +  1 

NCCJNiK'.M  =  1()"(K'+1)+1 
(  RlC.ll  1  llXil 

ri  Si  ll-  (Jlldi.l  U.NWIDI  I  I  III  N 
(  I  1  I  'l  ( Ol  l  l  lCll  N  I 

COIl'dC.N)  =  ("»Y 
N  =  N  +  1 

NCONdC.N)  =  l()’‘d(-l)+l 
(  RKill  1  COl  1  l  lCll  N  1 

('C)i;i([(',N)  =  CYI.R 
N  =  N  +  1 
Na)N(IC.N)  =  7541 
INDIl- 

(■ 

C  I  RON  1  AND  MACK  COl'FFIClIiN  l'S  DFI’HNDING  ON  WHICH  HDGE 
C'  FRONT  FDGF. 

IF'  (IHd).FOd)  THEN 
C  I  RON  1'  COEF'FICIEN  r 

COEF'dC.N)  =  CXFM 
N  =  N  +  1 
NCONdC.N)  =  7521 
(  MACK  COl  I  I  ICIFN  1 

COEFdC.N)  =  C\X 
N  =  N  +  1 

NCONdCNi  =  KEdC  +  NWIDl  )+l 
C  MA(  K  I  DGI 

I  1  SI  11  (111(1)1  O.NDI  I  Id  1  HI  N 
(  I  RON  1  COl  ITICIFN  I 

(  (d  Fd(  ,N)  =  (  XA 
N  =  N  +  1 

NCONdC.N)  =  lOdlC-NWlDl  )+l 
(  MACK  COE.FFICIEN  l 

COEl  (IC.N)  =  CXFM 
N  =  N  +  1 
NCONdC.N)  =  7531 
I  NDIF 
C 

C  rOF  COEFFICIENT 

COEF(IC.N)  =  E/,C 
N  =  N  +  1 

NCONdC.N)  =  10*(IC-NPE)+1 
C  MO  TFOM  COEFFICIENT 

COEFdC.N;  ^  CZB 
N  =  N  +  I 
NCONdC.N)  =  '’5(.l 
C 
( 

(  - ”***'1  RON  1  AND  MACK  EDtd-.S  1  XCEUDING  cokm  i^s'*********’**”******** 

C 

I  I  SI  IF  (dHd).i  (J.i.OR.IHd).FO.NDFMd-A.\D.(.IHdi.NF.!.OR 
+,IH(I).N1  .NWIDI  ))  Fill  N 
(  NL'.MMFR  OF  CONNl  (HONS  PER  NODE 
NCON(  IC.N  )  =  (1 
C 

(  El  Ed  (Ol  11  ICII  N  l 
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c:ol;f(1C',N)  =  cyy 

N  =  N  +  I 

NCON(IC.N)  =  10*(I('-1)+1 
(  KKill  r  ( OI-l'I  K'H-N  T 
COt-.I-dC.N)  =  CYY 
N  =  N  +  1 

NCONdC'.N)  =  1()*(IC+1)+1 

( 

(  I  RON  1  AM)  HAl  k  (  ()l  I  l  ICll  N  IS  1)1  I’l  NDINCi  ON  U.ilCll  FDGl; 

C  IKON  I  I  l)(,l' 

M  (111(1). l  U.l)  IIII  N 
(  IKON  I  COM  1  ICII  N  l 

COI  I'dC.N)  =  C.NIK 
N  =  .N  +  1 
NCONdC.N)  =  Cdl 
(  liACk  COI  ri  R'Il  N  l 

COM  (IC'.N)  =  C.N.X 
N  =  N  +  1 

NCONdC.N)  =  10*(lC  +  NWID[d+l 
C  RACK  MX,;-: 

M.SMl-'  (IH(I).f:a.NDFEP)  IHI-N 
C  I  RON  I  COM-'I-'ICIFNT 

COI  i-dc.N)  =  '::xx 
N  =  N  +  1 

NCONdC.N)  =  1()MC-NWIDH)+1 
c  ij.ACK  coF:r-’FiciE:N  r 

COM-dC.N)  =  CXFU 
N  =  N  +  1 
NCONdC.N)  =  CS.d 
i-NDIl- 
C 

C  lOldOM  I  ic  il-N  I 

COMilCNj  =  1  /(■ 

N  =  N  1 

NCONdC.N)  =  lOdK'-M'M+l 
(  HOriOM  COM  IK  II  NT 

COI-l-dC.N)  =  (  /li 
N  =  N  +  1 
NCONdC.N)  ■= 

C 

l-NDll- 

lo:  coNriNCT 

I  NDII- 
C 

(  tiL  NERAri;  DATA  FILE  VALUES 
(  TOTAL.  NODES  FOR  TdlS  ASPL.C  I  RA HO 
IL  (EAR  EO.0.0)  THL:N 
COUNT  =  .VNPL 
LISE 

COUNT  =  .KNPl.  +  :*NWIDE 
ENDIE 

( 

(  NUMBER  OI  C ONS  T.AN  T  TEMPI  R.ATURE  INPUTS 
CONH-MP  =  h 


90 


C  DUMMY  VAKlABLli 
ZER  =  0 
C 

C'  UNITS  TO  BE  USED 

IF  (UNITS.EO.  E  )  THEN 
USEE  =  1 
ELSE 
USEE  =  2 
ENDIF 
C 

C  PROBLEM  CAPABILITY  LINE 

c  maximum  nodes 

NMAX  =  750 

C  MAXIMUM  CONST.AN  r  TEMPERATURES 
I'MAX  =  50 

C  NUMBER  OF  lll.Al'FRS 
H  TRS  =  0 

C  DATA  SI  TS  Rl  OUIRED 
D1  =  2 
D2  =  4 
D3  =  0 
D4  =  0 
D5  =  0 
D6  =  0 
D7  =  0 
C 

C  ACCURACY  LINE 
C 

C  ACCURACY  BETWEEN  ITERATIONS 
ACC  =  0.05 

C  DAMPING  VALUE 
DAMP  =  0.66667 
C  MAXIMUM  ITERATIONS 
MAXIT  =  12 

C  CONVERGENCE  FACTOR 
CONFAC  =  0.8 
C 
C 
C 

C  CREATE  DA  FA  FILL 

OPEN  (. 4.FILI  .  =  N AMI, .FORM=  FORMA  FI  F:D',AC:Ci;SS=  DiREC  r-.RECL=  108 
+  .SLATUS=NFW) 

C 

('  LINE  1,  ri FLE 

WRITE(3,909)  DAFAJ- 
909  FORMAT!  IX.  A79) 

c; 

C  LINE  2,  PROBLEM  DATA 

WRITE(3.908)  COUNT, CONTEMP.ZER.ZER.ZER.ZER.ZER  ZER, USEE 
908  FORMAT(2X,9(I3.5X)) 

C 

C  ANALYZER  CONIROL  LINE 
WRITE(3,907)  ZER,ZER,ZER 
907  FORMAT(2X,3(I3,5X)) 

C 
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C  PROBLEM  CAP/\JilI.H  Y  L.lNl- 

WRrri;(3,908)  NMAX,  rM/YX,H  l'RS,Dl,D2.l)3,D4.D5.D6.i)7 

(' 

C-  ACCURACY  I.INI-: 

VVRIT{-(3.90.'S)  ACC.DAMP.MAXI  I.CONI  AC.I  I 
403  l-ORMAr(lX,2(l-9.7.1X),14X,i2.1X.l'4.7,lX,l-'4.5) 

C 

C 

C  CONS'UYN  r  rCMPCRA  rURf;  LINE 

WR1TE(3.906)  UPRT.LT.BLR'LLT.LWR'r 
900  FORMAT(lX,6(F12.9,lX)) 

C 

C  COEFFICIENT  EQUATIONS 

DO  112  1=1, COUNT  ^ 

WRITE(3,910)  (NCON(I,J),J=l,8) 

910  FORMAT(I4.3X,7(I4,8X)) 

WR1TE(3,911)  (COEF(I,N),N=l,7)  { 

911  FORMAT(7(F9.3,3X)) 

112  CONTINUE 

C 

C 

C 

CALI.  CLS 
WRriE(L994)  NAME 

994  FORMAK  //'.'  lilE  OU'ITUT  DATA  HAS  BEEN  PLACED  IN  A  I'lLE  NAMED  ' 

+,A6) 

END 


( 


r 
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